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X ■ We give results for the B and the D meson spectrum using NRQCD 

on the lattice in the quenched approximation. The masses of radially and 
orbitally excited states are calculated as well as S-wave hyperfine and P-wave 
fine structure. Radially excited P-states are observed for the first time. Radial 
and orbital excitation energies match well to experiment, as does the strange- 
non-strange S-wave splitting. We compare the light and heavy quark mass 
dependence of various splittings to experiment. Our B-results cover a range in 
lattice spacings of more than a factor of two. Our D-results are from a single 
lattice spacing and we compare them to numbers in the literature from finer 
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lattices using other methods. We see no significant dependence of physical 

results on the lattice spacing. 

PACS: 11.15.Ha 12.38.Gc 14.40.Lb 14.40.Nd 
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I. INTRODUCTION 



Mesonic bound states consisting of a single heavy quark, b or c, and a light quark, u, d 
or s, as well as gluons, provide an interesting laboratory to study strong interactions. The 
typical momentum within such states is much lower than the mass of the heavy quark. This 
leads to a situation where the heavy quark becomes non-relativistic and the properties of 
the bound state are essentially determined by the light quark and the glue. At leading order 
the splittings within the spectrum become independent of the properties of the heavy quark, 
such as its mass m.Q and spin sq, so that orbital and radial excitation energies are expected 
to match between the B-system and the D-system. The resulting approximate SU(2Nh) 
symmetry, with denoting the number of heavy flavours, is usually referred to as heavy 
quark symmetry, see and the references therein. At the next order, 1 /m.Q effects give rise 
to fine structure in the spectrum, several times larger in the D-system than for the B, see 
e.g. for a review. 

The spectrum of B and D states is not yet well established experimentally although 
several new results have been reported recently Here we study the spectrum theoret- 

ically and from first principles using lattice QCD. This will aid the experimental search for 
new states. In the case of well-established states it will provide a test for the theory and/or 
the systematic errors in our calculation. Of key interest are decay matrix elements for B- 
factory experiments. Knowing how well the spectrum has been obtained gives confidence 
that we understand how to simulate B and D mesons reliably. This is important for the 
analysis of systematic errors in matrix element determinations. 

To formulate heavy b and c quarks on the lattice, a naive discretisation is inappropriate 
since the lattice spacings currently available are not small compared to the Compton wave 
length of those quarks (m.Qa > 1). Presently there are two different formulations available 
to simulate heavy quarks, non-relativistic QCD (NRQCD) |p|,[10[] and the heavy Wilson 
approach For the b-quark on present lattices both approaches become essentially the 
same. However, in this regime, NRQCD is to be preferred since the inclusion of higher order 
correction terms is easily implemented. 

In this publication we report on our calculations of the B-meson spectrum for two different 



values of the lattice spacing a. Together with the results of [T2], which were obtained with 
the same methods at another value of the lattice spacing, we can investigate the dependence 
on a of our results. Physical results must be independent of a and hence we can perform a 
test of systematic errors inherent in our calculation. We find no such errors at a significant 
level. In addition, on our coarsest lattice, we were able to simulate the D-meson spectrum 
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and compare to results using heavy Wilson methods on finer lattices (where NRQCD does 
not work well since am c < 1 ). Early results on our coarse lattice have already been published 
131. 



m 



Section gives details of the simulations we performed and section [III] gives details of 



our fit procedure. Section |TV| gives our determination of the bare b and c quark masses. 
Section [V] discusses the behaviour of the splittings in the spectrum that we obtain. This 
includes fits to the dependence of the splittings on the mass of the heavy quark. Section |VJ 
compares the results in physical units at different values of the lattice spacing and with 
previous results as well as with experiment. Readers interested in our results for the physical 
meson spectrum could jump directly to this section. Section [VI]] contains our conclusions 
and our best estimate for the B-spectrum, based on the combined input from three different 
values of the lattice spacing. 



II. SIMULATION DETAILS 



A. Gauge field action 

Our calculation was performed on two sets of gauge field configurations, which were 
generated using the Wilson gauge action 

S G = 13 Y. [1 " 3 ReTr ( U x,-U x+ v^U+ + ^U+ J] . (1) 

X,|0.<"V 

This action has lattice artifacts of 0{a 2 ). For the bare gauge coupling (3, we used 5.7 and 
6.2. The lattice volumes and the number of configurations are given in table |I[ We will refer 
to these configurations by their respective (3-values. 



B. Light quark propagators 



The light quark propagators have been generated with the use of the Sheikholeslami- 



Wohlert action, also known as the clover action 114 



-qj;ic sw K Y IpxF-vp.xg-vpMj 



(2) 



On the configuration set with (3=5.7 the clover coefficient c sw is set to its tadpole-improved 
tree level value c sw = 1 .5667, as determined from the 4th root of the plaquette |l5j . This 
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reduces the lattice spacing artifacts in the light quark propagators to O{oc s a, a 2 ). At (3 = 6.2 
we used the non-perturbative determined value of c sw = 1 .6138, which removes the O[oc s a) 
artifacts from the light quark propagator as well [JT5] . 



In reference [jl?]] the light hadron spectrum at (3 = 6.2 has been calculated using the 
non-perturbative as well as the tadpole-improved tree level value for c sw - No significant 
differences in the meson and baryon spectrum could be resolved between the two values of 
c sw . From this we expect the difference between tadpole and non-perturbatively improved 
light quarks at (3 = 6.2 to be well covered by the size of the statistical errors in our case 
as well. This allows us to compare our (3 = 6.2 results to the tadpole-improved results at 



(3 = 5.7 and in reference [12 



For each value of (3 we used 3 different values for the hopping parameter k. The actual 
values are detailed in table y. The table also contains the values of k c and k s from the 



UKQCD collaboration fl7|]lql used in our calculation. The use of these values is appropriate 
for the analysis in terms of chiral extrapolations and scale setting that we have done. We 
also carefully include systematic errors from different chiral extrapolations and associated 
uncertainties in setting the scale. A recent re-analysis by UKQCD of their light hadron 



spectrum jy|] gives somewhat different values for k c and k s . Our errors encompass any 
changes this would produce in our physical results. 



C. Heavy quark propagators 

The typical momentum scale inside a heavy light meson such as a B or D meson is of the 
0(Aqcd)> which is small compared to the mass of the heavy quark. Therefore the mass of 
the heavy quark vtlq represents an irrelevant scale for the dynamics of the mesonic bound 
state and it is possible to simulate these states on lattices with a lattice spacing larger than 
the Compton wavelength of the heavy quark. 

In our simulation we use a non-relativistic expansion of the heavy quark Hamiltonian, 
which is known as NRQCD @JTg . 



H = H + 6H , (3a) 
D 2 

Ho := -- , (3b) 

2m Q 

5H := -c 4T ^-cr ■ B + c 27r ^(D E — E D) — c 37r ^cr ■ (D x E - E x D) 
2m.Q om.Q Sthq 

(D 2 ) 2 2 D< 4 ' (D 2 ) 2 , N 

- Cl ^ + C5a 24n^- C6a T6W5- (3C) 
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Please note that the rest mass term of H has been omitted, resulting in a shift of the 
Hamiltonian, which is discussed in section [TV]. In the case of a heavy-light meson the 
NRQCD expansion has to be organised in powers of Aqcd/itlq pPfl . Here this expansion 
is used up to ©((Aqcd/ttiq) )• We also include the p 4 term, which is believed to be the 
leading term in (9((Aqcd/tti.q) 3 )- The last two terms correct for discretisation errors from 
finite lattice spacing in respectively the spatial and temporal derivatives, n is a stability 
parameter used in the evolution equation (B). The matching coefficients Ci, . . . , Cg are set 
to their tadpole- improved tree level values fl5|j . 

With the Hamiltonian H and 6H the propagator of the heavy quark can be obtained 
from a Schrodinger-type evolution equation 

G t+1 = (l -a±6H) (l -a^H ) n U+(l -a^H ) n (l - a^SH) G t for t > 1 , (4a) 
d = (l -a±5H) (l -a^H ) n U+(l -a^H ) n (l -a^5H)4) x . (4b) 

With (b x we denote the source smearing function used on the initial time slice. At (3 = 5.7 
we use 20 different values for tuq in the range 0.6 < cutiq < 20.0 and at (3 = 6.2 we use 
10 values in the range 1.1 < am.Q < 6.0. Details, including the n values, are given in the 
tables [TIT] and [TV} For each value of |3 we performed 3 different runs. At (3 = 5.7 we label 
them A, C and S; for (3 = 6.2 they are labeled H, N and P. 

For the S-wave mesons at (3 = 5.7 we used up to three different smearing functions, 4>g,o> 
(J)g,i and $g,2, i n the different runs. These are convolutions of Gaussian functions for the 
light and the heavy quark with radii as detailed in table |V[ The configurations were fixed 
to Coulomb gauge. A local sink will be denoted with (J)l. In most cases our final (3 = 5.7 
results were obtained with both sink and source smearing. 

For (3 = 6.2 we use smearing for the heavy quark propagators only. In run H at (3 = 6.2 
we applied a hybrid procedure of Jacobi smearing |3T] and fuzzing [E2| . For runs N and P we 
fixed the configurations to Coulomb gauge. We used hydrogenic wave functions <bHg,i, 4^,2 
and 4>He,i for run N. The indices 'g' and 'e' denote wave functions of the ground and first 



excited state. The details are given in table [VI|. In the P run we used Gaussian smearing 
with two different radii, otq = 2.5 and 5.0. 

The spin operators applied to construct mesonic states with the correct quantum numbers 
are detailed in table 1 of reference 
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D. Lattice spacing 



In the quenched approximation one obtains different values for the lattice spacing, de- 
pending on the quantity it is determined from. This is expected to be caused by the strong 
coupling a s running differently in the real world and the quenched theory. 

We use the physical mass of the p-meson [|J to fix the lattice spacing. This procedure 
is justified from the typical gluon momentum in a B or D meson being of similar size to 
the momentum in a light meson such as the tc and p. Since heavyonium states probe a 
higher physical scale these are not appropriate to fix the scale for a heavy-light system in 
the quenched approximation. Using the p-scale should take care of most of the quenching 
effects. 

The determination of m p is complicated by the chiral extrapolation required, see reference 
24| for a review. At (3 = 5.7 we use the result of ||18||. The result of the linear extrapolation 



in the light quark mass m q is quoted as the central value and the deviation of the quadratic 
fit is treated as a systematic uncertainty. At (3 = 6.2 a linear extrapolation is reported 



in reference [17]. We treat the difference to the 3rd order extrapolation from p5l as a 



systematic uncertainty. The numbers are compiled in table |V11| . We use: 

(3=5.7: a" 1 = 1 .1 16(12)(^ 6 ) GeV, a = 0.176S(19)(±g 8 ) fm , (5a) 
(3 = 6.2: a- 1 =2.59(^ )(^) GeV, a = 0.0762(1?!) ( 26 ) fm. (5b) 

For comparison, table [VI I| also shows the lattice spacing as obtained from the string tension 0" 
and the bottomonium splitting Xt> — As a physical value for a we choose a result obtained 
from a potential model fit to the charmonium spectrum [ 25j| . The lattice numbers originate 



from f27]]28|| . These results are in agreement with the outcome of the m p analysis. As 
explained above, the bottomonium system probes a different scale and the values obtained 



using it do not agree with the result from light spectroscopy [29 



III. FITTING TECHNIQUES 
A. Parametrisations 

At (3 = 5.7 we used several different smearings at source and sink. For hadron correlators 
with a local sink, we applied simultaneous vector fits, requiring the fitted mass(es) to 
agree for all propagators: 
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(c|)L(t)l4)i(0)) = ^A i , 1 cexp(-mi t t), 1<i<m, (6) 



k=1 



A t , 1t =((|) L |i|> k )<^icl4)i). (7) 

In the case of sink and source smearing, we used simultaneous matrix fits. In matrix fits, 
the fitted amplitudes are constrained in their relationship with each other as well: 

n 

(*j(t)l4»t(0)> =^B] t ik B i)lc exp(-m k t), 1<j,i<m, (8) 

k=1 

Bi,k = (^kl*i) • (9) 



The fitting techniques are described in more detail in reference | 23|| . We found matrix fits 



to be more precise with respect to statistical errors. Due to the omission of the rest mass in 
eq. ([5]) the fitted mass is shifted with respect to the bound state mass. We denote the result 
of the fit as the simulation mass, m S i m . The determination of the shift will be discussed in 



section IV 



To extract mass splittings we applied two different procedures. One is to fit the masses 
as above, take their difference, and then calculate the error from the bootstrap or jackknife 
samples of the difference. With this procedure one can easily take advantage of using different 



smearings. In the case of a single smearing function, a ratio-fit provides an alternative |20 
For this one divides the bootstrap or jackknife samples of the two propagators and fits the 
outcome with an exponential ansatz. The mass shift cancels out of the difference in both 
procedures. 



B. Pseudo-scalar and vector meson 



On the (3=5.7 configurations the simulation masses for pseudo-scalar and vector mesons 
have been determined most accurately in run A. In this run we only used the smearing 
functions c()g,i and <&g,2- We found the double exponential matrix fit with sink and source 
smearing to deliver the most precise result. For the fit range we choose the initial time slice 
tmin two time slices larger than the first time slice delivering a reasonable x 2 - 111 general 
we choose the number of dropped time slices multiplied by the number of propagators used 
for the fit to be larger than or equal to the number of fit parameters. The reason for 
this procedure is as follows. The first reasonable value of x 2 is observed once the residual 
excitations are just masked by the statistical uncertainties, which allows for them to be 
still of similar size. Each excited data point can eat up one fit parameter. Dropping as 
many data points as fit parameters delivers a fit which is entirely dominated by statistical 
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fluctuations. The residual fit range dependence of those fits becomes negligible against the 
statistical uncertainties. We judge x 2 values resulting in a Q > 0.1 as reasonable, where Q 
denotes the probability of a fit having an even higher value of x 2 - The final result is given 
in table IVffll . 

In run H at (3 = 6.2 we only had one smearing function available. We extracted the final 
results from single exponential fits to the propagators with source and sink smearing. Their 
fit results turned out to be more precise than the ones from using a local sink. The final fit 
range was determined such that we observed a reasonable x 2 an d achieved stability of the 
fitted result against variation of the fit range. The results are displayed in table [IX|. 

In the run N we used hydrogenic wave functions of different radii. We generated smeared 
local and smeared smeared meson propagators. However no cross correlators, e.g. <J>Hg,i at 
sink and cJ) H e,i a ^ source, were calculated. Hence eq. (H) was inapplicable and we had to use 
vector fits in the case of smearing at sink and source as well. 

In double exponential vector fits to two smearing functions, we observed extremely low 
values of Q. We observed that this is connected to unfortunate statistical fluctuations on 
certain time slices. However the fit parameters turned out to be stable with respect to 
variations of the fit range. These fits will be discussed in subsection [VB| in more detail. To 
obtain a more precise result for the S-wave ground states, we resorted to single exponential 
fits to single propagators and compared the outcome for the different smearing functions. 
This is shown in figure [I] for the pseudo-scalar propagator at am,Q = 2.5. The octagons 
indicate the final result for each propagator, as determined from the Q-value after dropping 
two time slices. Within statistical errors all results are in reasonable agreement with each 



other. For the final result, which is also included in table [IX|, we choose the smeared- 
smeared <t>H g) i propagator. In the end these deliver the more accurate hyperfine splitting, 
due to superior noise cancellation between the pseudo-scalar and the vector meson state. 

In this context it is interesting to note that the propagators with local sink and <J>Hg,2 
source smearing plateau much later than the others, but the results are in agreement with 
those from other propagators. 

For the S wave states in the run P we only had the Gaussian smearing at the source 
with radius aro = 2.5 and local sink. Since these propagators plateau quite late, we used 
t m in = 16 for the final result, the error bars are not competitive with those above. Since 



they are needed for the later analysis of the P states we include them as well in table IX 



To describe physical bound states involving light u and d quarks, the results of tables [VT l I 
and FX have to be extrapolated in the light quark hopping parameter. On both sets of 
configurations, the difference between the critical and normal hopping parameter is smaller 
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than the uncertainty we assigned to k c in table |J |17j and we use k c in our extrapolations. 
The normal hopping parameter is the one for which the extrapolations deliver the physical 
m n /m p ratio. 

Due to the high statistical accuracy we achieved at (3 = 5.7 in the pseudo-scalar case, a 
linear ansatz in am q := — - ) in a full covariant fit to all three data points, results in a 
fit with x /d.o.f. > 8/1 for am.Q < 10. This corresponds to Q < 0.004. The resulting curves 
do not describe the data. We carefully checked whether this is caused by a residual fit range 
dependence and found all the fit parameters including the would-be strange to non-strange 
meson splitting to be stable against variation of the fit range. This was done for an initial 
time slice t min in a range from 3 to 6. 

We therefore extracted our final result from a linear spline to the points with highest 
and lowest m q and use the deviation of a quadratic spline as a systematic uncertainty of the 
chiral extrapolation. An example for the extrapolation is given in figure 0. From the figure 
it is obvious that interpolations to extract the heavy-strange meson mass are insensitive to 
the different ansatze and we do not assign an uncertainty due to the different interpolations. 
However, in the case of the heavy strange meson, we are faced with the problem that k s is 
highly sensitive to the quantity it is determined from. Our central value is interpolated to 
the k as determined from m^/m p , and the difference to the outcome for k corresponding 
to m^/rrip is treated as an uncertainty of the quenched approximation. The results are 
presented in table [X|. 

For (3 = 6.2 the statistical accuracy is not as high and our data are well described by 



linear extrapolations. The results are presented in table XI. 



IV. HEAVY QUARK MASSES 

A. Mass shift from dispersion relation 

The omission of the rest mass term itlq in the Hamiltonian eq. (^) causes most of the 
shift of the simulation mass m S i m with respect to the physical meson mass. The mass, 
m re i, of the meson can be determined from the relativistic dispersion relation of the meson 



E(p) = yTTlrel + P 2 ; which gives 

p 2 -[F(p)-F(0)F 
mrel= 2[E(p)-E(0)] • (10) 

Here E(p) denotes the total energy of the meson. The mass shift A re i is defined as the 
difference 
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A re i := m re i - m sim . (11) 

This shift per heavy quark should be universal for all hadronic states simulated at the bare 
heavy quark mass tuq. 

In our calculation at (3 = 5.7 we determined the mass shift from the difference in energy 
of the pseudo-scalar meson propagators with a|p| = and 2n/M. This was done in run 
C at k = 0.1400 with source smearing (bcj and a local sink. At large values of ttlq, we 
found a single exponential fit to the ratio of the correlators to plateau much later than the 
fits to the individual propagators. This is reflected in a large fit range dependence of the 
jackknife difference of the masses of the individual fits, for time slices in which no plateau 
was observed in the ratio-fit. For our final result we choose a minimal t- value two time 
slices larger than the first t-value for which we obtained a decent x 2 i n a fit to the ratio of 



propagators. The final result is presented in table [XII] and figure y. 

We also tried simultaneous vector fits according to eq. (|6]) with two exponents. We 
used propagators with source smearing (bcj and (J>g,2- The jackknifed difference of the 
fitted ground state mass is in agreement with the above procedure; however, the statistical 
uncertainties, especially for large values of ttiq, are larger. 

For (3 = 6.2 we calculated the mass shift in heavy quarkonia, since the statistical precision 
for heavy-light correlators at finite momentum was not sufficient. In the following, mass 



shifts from heavy quarkonia will be denoted by Ah- We simulated the vector-meson for 

'24" 



CL|p| < 2f£. The kinetic mass mi was obtained from fits to the dispersion relations: 



P 2 P 4 

Fsim(p) = m + - , (12a) 
Wt^mo + iL-^, (i 2b ) 
P 2 

Fsim(p) = TTLo + - , (12c) 

2m! 

with parameters mo, mi and m.2. E smi (p) denotes the simulation energy as determined from 
the propagator falloff. In the case of am.Q < 1 .3 we used the ansatze ( |12a|) and ( |12b|) ; for the 



three heavier tuq- values, ( |12b| ) and ( |12c| ). All fits gave fit parameters which were consistent 
within half of the statistical error. To obtain the shifts required for heavy-light spectroscopy 
we subtracted the simulation mass of the quarkonium vector-meson and divided by two. 
The final results are displayed in table [XIII] and figure |3|. It is interesting to compare to 



the result from reference - oAh = 1.29(2) obtained at arriQ = 1.22. Due to higher 

statistics, this result is much more precise. This value is included as a square into figure |3| 
and agrees well with the newer results. 
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B. Mass shift in perturbation theory 



The mass shift A can also be calculated in lattice perturbation theory [30 



A pcrt = Z m m Q - "E • (13) 

Here Z m denotes the renormalisation constant connecting the bare lattice mass m,Q with 
the pole mass and ho denotes the heavy quark self energy constant. In the perturbative 
expansion the 1-loop contributions from Z m and ho cancel each other to a large extent and 
the direct perturbative expansion of A pert is much better behaved than either perturbative 
series on its own. The Lepage Mackenzie scale aq* [I5| has been determined separately for 
A and it is larger than for Z m or for Eo. The coefficients for 

A pert = m Q [1 + <x s (aq*)-A (1) ] (14) 

can be found in table |XIV| . We use the (Xp(aq = 3.4) values as determined from the 1 x 1 
Wilson loop |31j with 2- loop running in order to evolve to the respective aq*. For the 
final mass shift we assign a relative uncertainty of ct^aq*). Since A' 1 ' is small, this is 
more conservative than the squared 1-loop contribution. The final results are displayed in 
table [XII| for |3 = 5.7 and table [XIII| for |3 = 6.2. For values of ttiq not included in table |XIV 
we interpolated linearly between the nearby values, which is completely sufficient within the 
claimed accuracy. The results for aA from perturbation theory and the lattice simulation 
are compared in figure 0. Apart from possibly the low ttiq region at |3 —5.7, the figure 
shows excellent agreement between the two ways of calculating the mass shift. 

For (3 —5.7 the stability parameter n differs in some cases between the perturbative 
results and the simulation. However for arriQ = 4, where perturbative results exist for 
n = 1 and 2, the effect of n is completely negligible: we obtain A pert = 3.88(22) vs 3.89(24). 
From a comparison of the simulation result of the runs A and C at (3 = 5.7 we can also obtain 
evidence of the effect of the different n on the simulation mass ra s ; m . For am,Q = 1 .0 and 
k = 0.1400 we measure am simiPS = 0.6265(21 ) and 0.6248(21 ) for n = 5 and 6 respectively. 
This difference is again completely negligible against the uncertainty we assign to aA, even 
if we enlarge it by a factor of 3 to allow for a larger effect between n = 4 and 5. The 
former n was used in the perturbation theory. Note also that this difference tends to be in 
the opposite direction to that in aA implying that the effect of n on the physical mass is 
reduced when compared to the shift. 

Here it is interesting to note that physical mass differences like the hyperfine splitting 
TTLhpf = m sim,v — TTi-sim.ps are even less sensitive to n. At the above mass parameter of 
am,Q = 1 .0 we measure am hp f = 0.0833(20) for n = 4 and 0.0835(20) for n = 5. 
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In summary the differences in n between the different runs as well as the perturbative 
shifts can be neglected safely even at the high level of accuracy we achieved here. This leaves 
us with a discrepancy between A pert and A re i for our lowest tuq- values, which is roughly twice 
as large as the uncertainty we assign to the perturbative result. 

On the other hand for (3 = 6.2 we observe excellent agreement between the precise result 
of with the perturbative calculation at the relatively low value am.Q = 1 .22. 



Given a value for the shift A and the simulation mass m S i m from tables IV11JI, 1X1 IXil 



and |Xnj, we can now calculate absolute masses for all the states. We do this for the ground 
state vector and pseudo-scalar mesons, both to fix the quark mass, as described in the next 
subsection. Moreover, we use the meson mass rather than the quark mass to discuss the tuq 
dependence, since it is more directly comparable to experiment. We frequently plot results 
against 1 /m,. aT , where Tru, av is the spin-average of the ground state vector and pseudo-scalar 
mesons 

m sav = J(3m v + m.p S ) . (15) 

This is preferable to m ps alone since the spin-averaging reduces the dependence on sub- 
leading spin-dependent terms. 



C. Bare heavy quark mass 

To determine the bare quark mass tuq corresponding to the b and c-quark, we compared 
the mass of the spin-averaged S-wave meson denoted with an overbar, with the experimental 
result. We used 5313 MeV for the B, 5405 MeV for the B s , 1973 MeV for the D and 2076 MeV 
for the D s ||. For the interpolations we used spline-fits to three neighbouring points. The 
fits were done quadratically in vtlq and 1 /tuq and no significant difference was observed 
between the two. From the strange and non-strange meson we obtained identical results for 
the quark masses: 

am b = 1.64(5)(l|Hi°), (3 = 6.2 (16a) 
am b = 4.20(25) (5) , (3=5.7 (16b) 
am c = 0.87[6) (3)(1? 3 ) , (3=5.7 (16c) 

The errors as indicated in the parentheses give the uncertainty arising from the mass shift 
and the statistical and systematic uncertainty of the a determination. The uncertainties 
associated with the simulation mass m S i m are completely negligible here. For m. b we used 
the perturbative shifts A pert . Using A H at (3 = 6.2 delivers am b = 1 -59 {^l 4 )^)^) and 
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using A rel at (3 = 5.7 gives cunt, = 4.1 6(1^) (7) (l^,), which is agreement with A pert but 
with larger error. For m c we used the simulation result A re i. Here A pert would give am c = 
1 .02(8) (2) (1° ). The deviation from the result eq. (|16c|) reflects the difference between the 
A-values at low arriQ discussed above. 

This careful analysis to fix the bare heavy quark mass is particularly necessary for fine 
structure splittings in the spectrum to be discussed in the next section. These are very sensi- 
tive to the heavy quark mass, generally as 1 /tuq. In addition, any errors in the heavy quark 
mass must be fed into errors in the fine structure splittings in order to avoid underestimating 
those errors. 

The bare masses do not scale with the lattice spacing as expected, because they are 
unphysical. A better quantity to consider would be the mass in the MS-scheme. This will 
be discussed in a future publication . 



V. MASS DEPENDENCE OF LEVEL SPLITTINGS 

In this section we describe how the results for the level splittings are extracted from our 
data. The dependence of the different splittings on the light and heavy quark mass is also 
discussed. 



A. Flavour dependent splittings 

The mass difference between heavy-light states distinguished only by their strangeness 
survives into the static limit. Based on the ideas of heavy quark symmetry such splittings are 
expected to depend weakly on the mass of the heavy quark. If a spin-averaged combination 
is taken, the leading heavy quark mass dependence arises purely from the kinetic term in 
eq. (|3b|). The size of the slope in 1 /m sav then gives information on the difference in (p£) for 
the b-quark in the strange and non-strange states. 

We calculated this splitting from the ground state S-wave results for k s and k c . The 
result is highly sensitive to the reported uncertainties in the chiral extrapolation and the 
determination of the strange hopping parameter k s . At (3 = 5.7 we determine the statistical 
uncertainties in a jackknife procedure applied to the difference of the individual results and 
at (3 = 6.2 we use the bootstrap. 

For (3 = 5.7 our statistical errors are very small and we consider additional systematic 
uncertainties for the chiral extrapolation and k s . For (3 = 6.2 the quality of our data is 
not as good and we give statistical errors only. Using the K s -value determined from the (£> 
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would lead to a 9% increase of the result, which is small compared to our statistical errors. 
The results are displayed in tables [XV| and [XV1| . In figure f| we plot the result for the 



spin- averaged splitting at |3 =5.7 versus the inverse of m sav , in order to display its heavy 
quark dependence. 

The figure displays a clear increase of the splitting with decreasing heavy quark mass. 
To quantify the slope of this dependence we perform a linear fit of the splitting result 
versus 1 /m sav . The result, converted into physical units, is detailed in table [XXV|. The 
slope corresponds to a (p£) difference of ~ [0.25(3) GeV] 2 , which is of the size of Aq CD , as 
expected. 

Because of the larger statistical uncertainties, we do not observe a significant slope at 
|3 = 6.2. The data can be fitted nicely to a constant. 

A comparison of our results with the ones of [|l^] for the pseudo-scalar case is plotted 
in figure ^. In this plot we show the result for the strange quark as determined from the 
K/p mass ratio only. Due to the large error bars we do not include the results obtained 
at the larger values of iuq for (3 = 6.2. Within the accuracy of around 12% in the case 
of (3 = 6.2 or better, no sign of scaling violations shows up in the plot. We also observe 
excellent agreement with the experimental result. 

B. Radial excitations 

In order to obtain a reasonably stable and long plateau for the radially excited S-states on 
the coarse lattice at (3 = 5.7 we applied triple exponential matrix fits to the three smearing 
functions 4>g,o> ^g.i an d 4>g,2- This was done for the run S for a single k of 0.1400 only, 
which is approximately equal to the strange as determined from the K-meson. Since in 
reference |L2| the dependence of the 2 1 So — 1 1 So splitting on the light quark mass was found 
to be very small, a variation of less than 2% when fixing k s from the K or K*-meson, we can 
ignore any mismatch in our k vs k s compared to the statistical uncertainties. We therefore 
treat our result as the answer for this splitting with k s as determined from the K. 

In figure ^| we show a typical example for the excellent stability of the simulation masses 
am sim against variation of the starting point t min of the fit range. The extent in t min for 
which we can resolve the excited state is 5 time slices or 0.28 GeV -1 . The rate of its 
disappearance is set by the 2S — 1 S splitting of 600 MeV. 

In the figure the first good value of Q is observed for t m ; n = 2. To be safe with respect 
to residual excitations we quote the final result for a fit range starting at time slice 4, which 
is the procedure described in subsection ITTT"Bl The peak in Q at t min = 5 results from the 
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fit becoming insensitive to the third exponential at this point. 

The results for all 6 heavy quark masses are given in table [XV1J| . In figure we plot 
the heavy quark mass dependence of the spin-averaged splitting. The result shows a clear 
increase as the heavy quark mass is reduced. In table [XXV| we detail fit results for the 
offset and slope of this splitting with respect to 1 /m sav . From the assumption that the 
increase of the splitting with 1 /m sav is caused by the difference in the kinetic energy p 2 /2rriQ 
between ground and radially excited states, the fitted slope gives a (p 2 ,) difference of ~ 
[0.95(15) GeV] 2 . This is of the size of a few times Aq CD as would be expected. 

On our fine lattice, since no cross-correlators between the different smearings had been 
calculated, we used simultaneous vector fits in all cases. The differences between the smear- 
ings <bHg,i, (t>He,i and 4>H g ,i turned out to be too small for simultaneous fits with three 
exponentials and we had to resort to fits using 2 exponentials. Again we choose the starting 
point t min of the fit range as described in subsection |IIIB| . With this procedure it is possible 
to extract reliable information on the excited state, as can be verified from the tables of 
reference [^3j for the case of T-spectroscopy. Using propagators with sink and source smear- 
ing in vector fits, leads to a suppression of the excited state contamination, which made 
it impossible to extract a signal for the excited state. Therefore we used propagators with 
smearing at the source and local sinks to extract the radially excited states. 

As denoted earlier, the fits to these propagators are plagued by statistical fluctuations, 
which lead to quite large x 2 and low Q values. However, the fits describe the data reasonably 
and the fitted parameters are stable against variation of t mm . This is shown in figure [5]. For 
a fit range starting point t m j n > 5 we obtain Q > 1 %, which is low compared to what we 
obtained in the other fits. However it is not that low, that the fit could be ruled out on 
statistical grounds. Together with the good stability of the fitted masses against variations 
of t m in as displayed on the right hand side of the figure, we believe that our signal is genuine. 
In this example we extract the final result from t m i n = 8. The results for the radially excited 



S-wave for this lattice spacing are compiled in table XVIII. 



C. Orbital excitations 



Orbitally excited P-state mesons have been investigated at both of our lattice spacings. 
The possible states consist of four non-degenerate energy levels; total angular momentum 
J = and 2 as well as two J = 1 states. In the heavy quark symmetry picture a jj coupled 
basis is appropriate. In the vicinity of the static limit, the J = 2 and the higher of the 
J = 1 states are close and separated from the J = and lower J = 1 . The former correspond 
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to a light quark total angular momentum of )i = I, the latter of jt = i. We use an LS 
coupled basis to study the states, but expect our 1 Pt and 3 Pi channels to mix, leading to 
the observation of the lower T p = 1 + state with both operators. We will denote the state 
corresponding to )\ = | with a prime. 

At (3 = 5.7 again we use one light hopping parameter, k = 0.1400. As in the case 
of the radial excitations, we treat this as the value corresponding to the strange quark as 
determined from the K-meson and the simulations have been performed in the run S. We 
used the derivatives of the smearings c|)g,o and (J)g,i at source and sink. The final results 
were obtained from double exponential simultaneous matrix fits to both smearings and are 



listed in table |XIX . 

The selection of the fit range proved to be very delicate for this a-value. The statistical 
error grows rapidly when increasing t m i n , since the signal to noise is exponentially related 
to the P — S splitting |5,33|. We give an example in section |V b] , where the fine structure is 
discussed. For the lightest values of am.Q and correlators with 3 P 2 wave operators in the T- 
representation we obtained small values of Q of a few permille. We include the corresponding 
mass values into the table for the sake of completeness and mark them in italics. However 
we disregard them in the further evaluation. The results are always in agreement with the 
ones obtained in the E-representation and we do not believe there is a serious problem with 
this, simply statistical fluctuations. 



In the second to last column of table [XIX| we give the spin-averaged P-state result which 
we calculate as 



m(P sav ) = — 



1 •m( 3 P )+3-m( 1 P 1 )+3-m( 3 P 1 )+5-m( 3 P 2 E)l . (17) 



The result is also shown in figure || For comparison we include the experimental result for 
the B*j(5850) resonance and the spin-average of the D s i and D* 2 0. The figure displays at 
most a mild heavy quark mass dependence. To quantify this, we report in table [XX V| on 
the offset and slope of this splitting in physical units. The slope is consistent with a (p£) 
difference of 0[Aq CD ), but is also consistent with zero. 

P-states were also investigated in the run P for (3 = 6.2. We choose am,Q = 1 .6, directly 
corresponding to the b-quark in eq. (|16a ). For the light quarks we use the strange hopping 
parameter k = 0.1346. Since the results on our coarse lattice depend only weakly on ttiq 
we take the outcome as the final answer for B s . 

In the simulations we used two different smeared sources together with local sinks. Again 
we used derivatives of Gaussian smearing functions. The masses were extracted from double 
exponential vector-fits to both propagators. We observe reasonable Q-values for all applied 
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operators and include all channels into the spin-average 



m(P sav ) = — [l • m( 3 P ) + 3 • mC?i ) + 3 • m^P, ) + 2 • m( 3 P 2 E) + 3 • m( 3 P 2 T) 



(18) 



The results for the fitted masses are displayed in table [XX| . The splitting between the 
spin-averaged P and S waves is given in tables [XX]| and |XXI1| . 



D. Radially excited P-states 



Having available 2 different smearing functions at both of our lattice spacings, it is 



possible to obtain information on the radially excited P-states as well. In figures [TC] we show 
the dependence of the fitted masses of the spin-averages of the 1 P and the 2P on the starting 
point t min of the fit range for two different values of tuq. The figure displays a clear signal 
for an excited state and reasonable stability with respect to variations of t m i n . However 
the error grows rapidly with increasing t mm . The Q-values of the 3 P2E-fit, which is the 
last of the individual states included in the spin-average to reach a plateau, are included in 
figure pj]. Decent Q-values are observed for t rmn = 3. Since this is a 6 parameter matrix fit 
to four propagators, we take our final result from t min = 5. The results for the spin-averaged 
2P-state are summarised in table [XIX]. We give the splitting to the spin-averaged IS and 



1 P-states in table |XXJ. We do not observe a significant slope for the splitting with respect 
to l/m sav . 

For (3 = 6.2 we show the plateau in figure 11. Due to the finer lattice, the growth in 
the error with increasing t m i n is smaller than before. For the mass parameters used here, an 



example of a Q-plot will be given in figure [TBI below. Here the first decent Q is observed for 
t m i n = 2. Since this is a vector fit to 2 propagators, we take our final results from t min = 6. 



The result and splittings are included in tables XX and XXII 



As noted when discussing radially excited S-states, due to our conservative selection of 
the fit range, we expect residual excitations to be negligible within the quoted statistical 
errors. 



E. Hyperfine splittings 

The mass difference between a pseudo-scalar and a vector S-wave meson is caused by 
the spin of the heavy quark. This hyperfine splitting is expected to vanish in the limit of 
infinitely heavy quark mass. 
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On our 6 = 5.7 lattice we determined the hyperfine splitting m.h p f from the difference 
of the results in table |V111| . A crucial ingredient in obtaining a small statistical error is to 
choose identical fitting ranges for both correlators. If one of them has a plateau at a larger 
value of t m in than the other, we took this larger value to obtain ruh p f from the difference of 
the fitted masses. The error in this procedure is estimated with a jackknife. The results are 
displayed in table [XX111| . 

The chiral extrapolation of the hyperfine splitting turned out to be less difficult than that 
for the simulation mass of the pseudo-scalar and vector mesons. The curvature seems to 
cancel out between them and we have been able to perform linear fits to extrapolate to the 
chiral limit. However, in order to be consistent, we assign a systematic uncertainty to the 
result. This uncertainty is obtained from the difference to the outcome of first extrapolating 
the individual mesons to the chiral limit and then determining the hyperfine splitting. In 
this case we use the quadratic extrapolations to the chiral limit to take a possible curvature 
into account. 

In our chiral extrapolations of the hyperfine splitting we observe a negative slope with 
respect to the mass am.q of the light quark, which is illustrated in figure |T2|. The left 
hand side gives an example of our chiral extrapolations and the right hand side shows 
the slopes measured at each value of am.Q. In order to construct a physically meaningful 
quantity, the latter has been multiplied by the strange quark mass, such that we can compare 
with experimental results for the difference between the strange and non-strange hyperfine 
splittings. 

For B-mesons the light quark dependence of the hyperfine splitting is not well resolved 
experimentally, because of large uncertainties in the B s hyperfine splitting. For D-mesons the 
situation is much clearer and one observes an increase with the light quark mass. However 
the magnitude of the slope is largely dependent on whether you compare the D s hyperfine 
splitting with the hyperfine splitting of the D + or the D°. We expect this difference in the 
experimental results to be mainly due to QED-effects, since these come in with opposite 
signs in the D + and the D°. Since the D s is positively charged as well, QED effects should 
largely cancel when comparing the hyperfine splittings of the D s and the D + and one obtains 
a positive slope for the D-meson from the experiment. 

Comparing our data to the experimental results, one observes our hyperfine splittings to 



be too small. This will be discussed in more detail in section |\T|. With respect to the slope, 
the result at the D has clearly the wrong sign and its magnitude is approximately twice as 
large as that from the experiment. We did not observe this effect in our 6 = 6.2 results, 



neither was it observed in [|12 |. Both of these results did not achieve the high statistical 
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accuracy we have at (3 =5.7 and also use values of the heavy quark mass at around the B 
or heavier. For those values of trq the light quark mass dependence at (3 = 5.7 is also not 
that significant. 

A slope of similar sign and size has been observed in the calculations of [|34],BH] , although 



the authors did not comment on this. Reference 34] used a highly improved gluonic action 



with NRQCD heavy quarks on even coarser lattices and a heavy clover action for the 
heavy quarks on a finer lattice (3 = 6.0. A detailed comparison with these results will be 



given below in section [VT|. In this context it is interesting to note, the slope of the hyperfine 
splitting as a function of the quark mass turns out to be to small in light hadron spectroscopy 



as well 44 



The calculations listed above are performed in the quenched approximation, which could 
be a factor contributing to the wrong slope. In potential model language, which is not 
necessarily appropriate here, the hyperfine splitting is related to the square of the wave- 
function at the origin. This in turn depends on the light quark mass and is independent of 
m.Q as rriQ — > oo. The wrong slope could then reflect the fact that the wave-function at the 
origin is not increasing rapidly enough as the light quark mass increases. This is natural 
in the quenched approximation as the potential at the origin is weakened by the coupling 
constant running too quickly to zero at short distance. 

An alternative scenario is one in which the coefficients of the relevant terms in the action, 
here c sw in the light quark action, effectively carry some quark mass dependence that has 
not been included, leading to an underestimate of the hyperfine splitting at large am q . In 
this case the effect would disappear as a is reduced and this seems to be contradicted by 
results on finer lattices . 



Another cause could be a problem in the chiral extrapolation itself. The experimental 
result for hyperfine splitting J/ih— r| c in the charmonium system is smaller than the hyperfine 
splitting for the D and D s -mesons. If one considers charmonium as a D c -meson, one has 
to conclude that there is a maximum of the hyperfine splitting as a function of the light 
quark mass for m q < m c . If this maximum is attained for m q < m s , our observation of a 
negative slope of the hyperfine for m q ~ m s would be in agreement with nature. In this 
case extrapolations from the strange region to lighter u and d-quarks as well as the chiral 
hmit would be impossible. 

Our final results for the hyperfine splitting at k c and k s are given in table [XXIII] . In 



figure [13] we display the dependence of the hyperfine splitting on the spin-averaged heavy- 
light meson mass. A linear fit in vci^ for the five heaviest m sav values gives reasonable 



values of Q. In table |XXV| we give the numerical outcome of this fit for the strange and 
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non-strange hyperfine splittings. As expected from HQET the intercept always turns out to 
be zero within statistical errors. 

In the H run for |3 = 6.2 we determine the hyperfine splitting from ratio-fits. In order to 
determine this without excited state contamination we use a fit interval for which both of 
the individual correlators have reached a plateau. As noted above, no significant dependence 
on the light quark mass was observed and we were able to fit the results to a constant with 
reasonable values of Q. The result is given in table |XX1V| . 

In case of the N and P run we determined the hyperfine splitting from the jackknife 
difference of masses obtained from the pseudo-scalar and vector meson propagator. In the N 
run we compared the outcome for the different smearings available, for different values of the 
starting point t min of the fit range. An example is shown in figure [TR The different smearing 



functions lead to compatible answers for the hyperfine splitting. We use the outcome from 



the propagators with sink and source smearing 4>Hg,i for our final result. In figure [L6| we 
compare the outcome of the different runs at (3 — 6.2. Clearly the outcome from the run N 
is the most precise. The result for the physical B s hyperfine splitting will be extracted from 
this results. 

Having observed clear signals for the radially excited S-wave states on our coarse lattice, 
we also studied their hyperfine splittings. Unfortunately the statistical noise grows rapidly 



and we observe no clear signal for a non-zero splitting. Our results are given in figure TA, 
comparing the radially excited state hyperfine splitting to that of the ground state. Although 
we cannot give a value for the hyperfine splitting of the radially excited S-state, our results 
support the expectation that it should be equal to or smaller than the ground state splitting. 



F. P-state fine structure 

To extract the P-state fine structure we investigate the jackknife difference of the masses 
of the individual channels reported in tables |X1X| and [XX|. Because the statistical noise 



grows rapidly as t min increased, this proved to be delicate. For (3 = 5.7 this is illustrated in 
figure |17| For the matrix-fit to the 3 P2p propagators we observe a jump in Q for t m j n = 3. 
However in the plot of the fit range dependence of the 3 P2p — 3 Po splitting the statistical 
uncertainty doubles between t min = 3 and 5. Therefore it is hard to tell whether there is a 
plateau or not. 

We quote final results for t m j n = 5, which corresponds to dropping 2 timeslices from 
the first reasonable Q-value. With this procedure we obtain a large statistical error and 
no significant splitting can be resolved. More aggressive fitting would have led to a result 
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compatible with zero but with a statistical error of ~ 30 MeV. We give our final numbers 
in table |XX1| . For the splitting we used the same fit range for both channels, which leads to 
slight deviations from the direct difference of the results in table £XTX. 



For (3 = 6.2 the situation is easier, as shown in figure |18[ The noise on the splitting does 
not grow as fast as on the coarse lattice, because the P — S splitting is smaller in lattice 
units. We observe the first reasonable Q values for t min = 3. Since this is a 6 parameter fit 
to two propagators, we drop 3 time slices and quote the final result for t min = 6. The results 
are given in table |XXII| . Here we also quote results for the splitting of the J = 1 channels 



to the J = state. We reiterate that no significance should be attached to any difference in 
the results between the 3 Pi and the operators. 

VI. THE PHYSICAL MESON SPECTRUM 

In this section we determine the physical B and D-meson spectrum and investigate 
scaling by comparing results at different values of the lattice spacing. We also compare with 
experimental results and other lattice calculations. 

A. B-meson spectrum 

At both of our lattice spacings we can simulate the b-quark directly. Here we discuss 



our results for the physical B-spectrum. Together with the findings of [12] we want to 
investigate the dependence of the individual splittings on the lattice spacing. The findings 
are compared both to the existing experimental results and lattice investigations performed 
by other groups within a similar framework using NRQCD |53].|33 . 



The results of reference [|34| are useful in that they work at a larger lattice spacing than 
we do here. There are a number of problems, however that make their results not directly 
comparable. For example, they do not use either smeared correlators or standard fitting 
techniques, and this will give rise to an unknown systematic error. In addition they do not 
see a difference between fixing the lattice spacing from m p using the clover action and from 



the charmonium 1 P — 1 S splitting. This is clearly seen on finer lattices [[L7],|1^,|3(J . If this 
arises from overestimating a -1 from m p because of discretisation errors, then this is another 
source of systematic error. In particular, this feeds into the fixing of the bare b or c quark 
mass and into hyperfine splittings. Their final result for the splitting does not take into 
account the effect of any of the uncertainties in the bare quark mass determination. This 
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is particularly important for the hyperfine splitting and causes their errors to be heavily 
underestimated. 

The results of overlap with, but are not as complete as ours. 

Unfortunately there are no results for heavy clover fermions available that we can use. 
References |]38| , |39[| quote numbers for the B spectrum. The first reference still uses extrap- 
olations from the lighter quark masses into the b-region. Reference |!9| determines the 
bare b-quark mass from heavyonium, which is not suitable for the heavy clover approach 
at the lattice spacings used f|0],f|l]. However we will later compare to their findings for the 
D-spectrum, since this problem is not so severe for charmonium at the lattice spacings used. 

Results from taking the b-quark as a static source also exist for spin-independent and 
flavour splittings, which survive in the infinite mass limit, see for example JIIJ. However we 
restrict the discussion here to a comparison with results simulated at the physical b-quark 
mass directly. 

In the following we denote spin-averaged states by an overline. 

We summarise our results for the B-spectrum in the tables [XXV1| and [X.XV11| . As an 
example of the splitting between a strange and a non-strange meson we discuss the difference 
of the pseudo-scalar B s — B splitting in figure We observe no scaling violations between 
our results and the agreement with the experimental value is excellent. 



The results from [35] given in figure [T^ include the statistical errors only and are taken 
from their figure 15. In [3~5| additional uncertainties for the average of the results at the two 
finest lattices are mentioned in the text. The overall agreement with our results is good. 
They notice an upward jump, however, for their result on their finest lattice. Unfortunately 
our result from our finest lattice comes with large uncertainties so that we are unable to 
clarify whether there is any real effect here. Given the lack of scaling violations in the rest 
of the results, it seems unlikely to us. Table |V11| confirms that the scale from m p used by 
us and the scale from c used in |35| are very close. The mismatch of scales of ~ 3% can be 
neglected safely. 



The results from [M\ are also in agreement with ours. They use the K*/K ratio to fix 
the strange quark mass. This reduces their results compared to that using the K/p ratio. 
Assuming a shift of 10 to 20 MeV from this would increase the agreement. This is the size 
of the effect we observe on our coarse lattice from fixing the strange from cb/p instead of 
the K. 

Figure |2(] shows the scaling of the radially excited B s -meson. As discussed in subsec- 
tion |V Ij| already, the extraction of a result for our finest value of a turned out to be more 
problematic than anticipated, and we are left with quite large statistical uncertainties. How- 
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ever the final result is in good agreement with the result from our coarser lattice as well as 
the result of JIJ. 

We also included a preliminary experimental result from the DELPHI collaboration for 
an admixture of the non-strange B' — B and B*' — B* splitting f|],f|. Assuming the hyperfine 
splitting of the two states to be of similar size, which our findings support, we observe 
reasonable agreement here. Table [XXVJ| also contains results for the radial excitation energy 
of the vector state and the spin-averaged S-wave. 



The orbital excitations are compared in figure [21]. The lattice results for the spin-averaged 
strange P state scale very well. The magnitude agrees nicely with the B*j(5850) resonance, 
which is expected to be an admixture of the two )\ = I states. 



Our results for the radially excited P-states are compared in figure g^. This is the first 
observation of a signal for these states in a lattice calculation. As the figure shows we get 
consistent results from the two different lattice spacings investigated. To the best of our 
knowledge radially excited P-states have not been observed yet experimentally. 

The splittings discussed above are all essentially light quark quantities, which survive 
into the static limit. Their scaling or non-scaling says more about the light quark action 
than the heavy quark sector. The hyperfine splitting is of a different nature and from its 
scaling behaviour one can learn about how well the heavy quarks are being described on the 



lattice. We display this in figure |23[. Our result for the strange hyperfine splitting together 
with the findings from |T2] shows good scaling. 

However the result is much smaller than the experimental value. Since the leading term 
in the hyperfine splitting arises from the a ■ B term in the action, eq. fl3c|) , the result for the 
splitting is sensitive to the coefficient C4 and the inclusion of radiative corrections beyond 
tadpole improvement is required. Preliminary calculations indicate that the inclusion 
of the 1-loop corrections would increase the hyperfine splittings on the order of 10% for 
the lattice spacings used. The quenched approximation might also play a role here, since 



in light spectroscopy the hyperfine splittings turn out to be too small as well, see HJ] for 
a review. This effect increases with increasing quark mass. Unfortunately reference p? 



which investigates the effect of the inclusion of two flavours of dynamical quarks on the 
B-meson, does not give any evidence for an increase of the B* — B splitting due to sea-quark 
vacuum polarisation effects. 

From the experience in T-spectroscopy using NRQCD, one could have expected to 
observe scaling violations in the hyperfine splitting. Using the same heavy quark Hamilto- 
nian as we do, ]2D|] reports an increase of 50% for the T — r|t, splitting, within the range 
from (3 = 5.7 up to (3 = 6.2. The leading discretisation correction for the hyperfine splitting 
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is C((apgi U on) 2 ) [jnj- Typical gluon momenta for the T-system are ~ 1 GeV, while for the 
B-system they are C(Aq CD ). From this one expects reduced scaling violations of ~ 10% in 
the B-system for our range of lattice spacings. This is the same size as our uncertainties on 
the hyperfme splitting and therefore consistent with the fact that no scaling violations show 
up in figure p3l 



Results from |35| are for the chirally extrapolated splitting B* — B with statistical errors 
only, taken from their figure 17. In the text, the authors quote a result for the strange 
hyperfme splitting from the average of the two finest lattices, which is 3 MeV higher than the 
same average for the non-strange hyperfine splitting. An upwards shift of 3 MeV increases 
the already excellent agreement even further. 

The results of ||34||, on the other hand, exhibit a clear disagreement to our findings as 



well as the findings of ||35|| . We believe that this is because they have determined the bare 
b-quark mass arrib using heavyonium. 

The fine structure of the P-states is the last topic of this section. Unfortunately we have 
not been able to resolve this clearly on our coarsest lattice. The situation is displayed in 
figure 123, for the three sublevels which were resolved at (3 = 6.0 and 6.2. To investigate 
whether there is evidence for scaling violations in the fine structure we calculate the jack- 



knifed difference of the highest and lowest state. This is shown in figure |25|. The error bars 
turn out to be large and the figure is inconclusive. A more aggressive fit on the coarser 
lattice, as discussed in subsection |V E| would lead to the conclusion that scaling violations 
were seen, but we believe that further work is needed to resolve this question. 



Our results for the B-meson spectrum, together with those of [12j do not show signs of 
residual lattice spacing dependence within the achieved accuracy. Therefore we can average 
the results for the different values of the lattice spacing a to obtain our final results on 
the quenched B-meson excitation spectrum. The averages were obtained in the following 
way, for each value of a we add the different uncertainties in quadrature to obtain a single 
value. Here we omitted those sources of uncertainty which are associated with the quenched 
approximation. These are the uncertainty arising from fixing the strange quark mass from 
different physical quantities and in the case of the results of |12j the additional uncertainty 
of the lattice spacing a associated with the physical quantity used to fix a. At this step 
we also symmetrised with respect to unsymmetric uncertainties. The central values have 
been obtained from fitting the results to a constant in an uncorrelated fit with the above 
described uncertainties. This puts more weight on the more precise results than a simple 
average. 

Our analysis at the individual values of a does not include an uncertainty for the residual 
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effect of the lattice spacing. Our actions are improved to O{oc s a, a 2 ). Therefore for each 
value of a we add the maximum of a s aAq C D an d cl 2 Aq CD in quadrature to the uncertainty 
used in the fit. We quote the smallest of these three so obtained uncertainties as our final 
uncertainty for the quenched B-meson spectrum. This way we quote an accuracy which is of 
the same size as the one we checked for scaling violations. Determining the final uncertainty 
from the x of the fit would reduce the uncertainty beyond this level. This procedure also 
ensures that residual lattice spacing artifacts are properly included if the achieved accuracy 
differs over the three individual results and the average is largely determined by the coarser 
lattices. 

Our final result on the B-meson splitting spectrum in the quenched approximation is 
given in table |XXV111| and figure EBI. 



The question of the effect of quenching on the spectrum goes beyond the scope of this 
paper. We refer the reader to |}7[] . There the effects of the inclusion of 2 flavours of dynamical 
quarks on the spectrum in NRQCD have been investigated and compared to the findings of 
||12|| . No significant difference between the quenched and rif = 2 results for the IS, 2S and 
1 P-states was found. In particular, as mentioned earlier, no significant sea quark effects were 
seen in the 1 S-hyperfine splitting. Since our investigations confirm scaling in the quenched 



heavy-light spectrum the conclusions of [37fl are unchanged. 



B. D-meson spectrum 

In this section we discuss the D-meson spectrum and compare our result to existing 
lattice results as well as to experiment. 

The convergence of the NRQCD expansion is particularly important in the D-range, 
where the expansion parameter Aqcd/itlq ~ \- Useful results on the question of the 



convergence are contained in |34|. There the authors study the complete NRQCD ac- 
tion to C((Aqcd/itlq) 3 )- Here we include all terms up to (9((Aqcd/ttiq) 2 ) and the rel- 
ativists correction to the kinetic energy in O^Aqcd/ttiq) 3 ). The authors of [34| calculate 



the heavy light kinetic masses using eq. (|12dJ and show that the difference arising from 
C((AQCD/m.Q) 3 )-terms i s consistent with the expectation that they are sub-sub- leading in a 
Aqcd/ttt-q expansion. The changes to the spin-averaged meson mass that we use to fix the 

4 

quark mass are dominated by the -t- relativistic correction that we include. From this we 
conclude that remaining (^((Aqcd/tti-q) 3 ) and higher order terms in the NRQCD expansion 
would only change the physical masses by at most a few percent. This allows us to use the 
results of p| to estimate the changes in the hyperfine splitting which would be produced 
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by these additional terms at fixed bare quark mass. 

The authors of [Q find that the 0((AQ CD /rriQ) 2 ) terms produce an effect somewhat 
smaller than a Aqcd/tuq expansion might suggest, since they affect the hyperfine splitting 
indirectly. The only spin- dependent term at C((Aqcd/tuq) 2 ) is a spin-orbit type interaction. 
At C((Aqcd/iti-q) 3 ) most terms produce a change of a few percent, but the one which is 
directly related to the <r- B-term: {D 2 , cr- B} reduces the hyperfine splitting at the charm by 
20%. However when including the other operators of 0{{Aq CD /vtlq) s ), the second largest 
effect comes from the cr(ExE + BxB) operator, which is also spin-dependent and works in 
the opposite direction to the other one. The total effect of (9((Aqcd/tiq) 3 ) is below 10%. 
This is of the size of the naive expectation for the suppression with respect to the leading 
term and not at all inconsistent with good convergence of the NRQCD expansion. Since 
we do not include these terms, we conclude, that we may be overestimating the quenched 
lattice hyperfine splitting of the D-meson by 10%. 



Our results on the D-meson spectrum are summarised in table |XXIX| and in figure |27. 
The overall agreement to the experimentally observed spectrum is good. We will now dis- 
cuss the individual splittings in more detail. We also compare our results to the lattice 



studies of p^ , |38|j39|| . All of these results use the quenched approximation as well. The pub- 
lications D38],H| apply the heavy clover approach |11|| , which has quite different systematic 
uncertainties from NRQCD for charm quarks. 

The flavour dependent D^*' — D'*' splittings are in good agreement with the experimental 
results. Here it is interesting to note that our results reflect the increase of ~ 10 MeV from 
the B to the D-meson system, which can already be expected from the good agreement of 
the slope with the experimental outcome in table |XXV| . 



In figure |28] we compare our result for the strange to non-strange spin-averaged splitting 
to other lattice calculations. In order not to disguise other possible effects, we excluded the 
uncertainty of the strange quark mass from the plot. The results from |34j] use the K*/K ratio 
to define the strange quark. This should shift the results downwards, compared to fixing 
k s from the K/p ratio as used for the other results. The implications have already been 
discussed in the previous sub-section |V1 A| . The uncertainties again allow for an upwards 
shift of these results by 10 to 20 MeV. It should be noted that we combined the results 



for the hyperfine splittings and the pseudo-scalar D s — D splitting from [34| to obtain the 
spin-averaged splitting. 

Because of different systematic uncertainties, it is particularly interesting to compare 
to the heavy clover results of [PS[ . The results obtained with the use of the m p -scale, 
which is the same as what we use, agree very well with ours. This agreement is expected 
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since 



[38H uses the same light quark and gauge field action and this quantity is essentially 



determined by the light quarks and the gluon field. We conclude that the results for the 
flavour dependent D s — D splitting agree well between the different approaches and, within 
the accuracy achieved, agree well with the experimental result. 

For radially excited D^*' '-mesons no experimental results are known to us. However the 
DELPHI collaboration reports on the observation of the non-strange D*' |J. This result is 
still awaiting confirmation by the OPAL and the CLEO collaboration and its interpretation 



is disputed on the ground of its small experimentally observed width [p3| , f46f . The splitting 
between the DELPHI result and the D* has a similar size to our D*' — D* splitting. Ref- 
39| reports lattice results from the heavy-clover approach. From their plot we read 



erence 



D' s — D s « 840(160) MeV, which is in agreement with our findings. However this includes 
what the authors call " continuum" extrapolation out of a regime where the expansion pa- 
rameter am,Q = 0[\) is not small. We would prefer to compare to the unextrapolated 
results at the individual values of a. 

Experimentally the only well established charmed P-states in the particle data book || 
are those which are expected to correspond to the states of total light angular momentum 
ji = |. Recently the CLEO collaboration []7||| claimed the observation of the Di state 
corresponding to )\ = j. CLEO gives a preliminary result of D T = 2461 (+34) ( 1 0) (32) MeV, 
which is slightly heavier but compatible in error bars to the = 2425 MeV 0Q Our 
lattice calculation delivers the mass of the lighter of the two states. We did not observe a 
signal for an excited state slightly heavier than this. 



In table [XXIX| our result for the D s ( 1 P) — D s splitting is compared to the spin-average of 
the and D* 2 , the ji = § states. The agreement is reasonable. Reference |58] reports on 
the D s i — D s splitting from a lattice study with the heavy clover approach. A comparison to 



our result for this splitting is given in figure [29|. When using the same scale obtained from 
m p both lattice results agree very well with each other. The agreement with experiment is 
also good. 

It should be noted however, that the experimental result included in figure ^9| is not 
necessarily the same as ours. The experimental P-state corresponds to )\ = I. If the CLEO 
trend is confirmed and the Di is indeed heavier than the and the same holds for the 
D sl states, then the lattice result also corresponds to j\ = |. If not the lattice result will 
correspond to j\ = \, but the two states will be so close, that any mismatch is well covered 



We quote the charge-average. 
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by the error bars. 

Table [XX1X| contains our final result for the radially excited 2P-state. This is the first 
result for this state from a lattice simulation. 

As in the B-system the hyperfine splittings are too small when compared to the exper- 
imental result. Whereas in the B-system they were too small by ~ 40%, here they are 
low by « 25%. This could reflect a more severe quenching error for B-mesons. B-mesons 
are somewhat smaller states than D mesons and probe slightly different scales. This is a 
sub-leading effect in a heavy quark symmetry picture, however. Alternatively, if the error 
comes from radiative corrections to the C4 coefficient, that would need to increase with thq. 



That is seen by the authors of [43 



The large uncertainty of ~ 20 MeV on our result arises from the chiral extrapolations 
used in the lattice spacing determination and the way in which this feeds into the fixing of 
the bare quark mass. Naively we expect a doubling of the relative error, because a larger 
value of a requires a smaller value of am.Q to deliver the same physical m sav . This smaller 
value of qiuq gives a larger hyperfine splitting am,h. P f. When converting to physical units 
it picks up the uncertainty of a for the second time. In fact a factor of four is seen because 
of the flattening of the relation between the mass shift A and bare heavy quark mass thq, 
see figure 3, as well as the steepening up of the hyperfine splitting curve for large values of 



rrisav in figure [13. 

In figure [31] we compare our results to the results from |J3| obtained in NRQCD. We 
choose their result in 0((Aqcd/tti-q) 2 ) as most relevant for this comparison. The good 
agreement with our results is in fact misleading. They fix their c-quark mass from charmo- 
nium instead of the D. Fixing from the D would lead to a larger value of am c and lower 
hyperfine splitting, see subsection |1V C| and reference . 

It is interesting to compare the result for the hyperfine splitting between NRQCD and 
heavy clover quarks. This is also done in figure [30] for the D s hyperfine splitting. The 



heavy clover results of |3J| appear to be higher than our result. However this is a result of 
their higher choice of scale coming from J/i[> instead of from m p . This is confirmed by the 



findings of Using a scale from m p gives a result which agrees with ours, using a scale 
from J/il) agrees with reference [B3. It should be noted that in [3~8| the bare quark mass is 
determined from the D, where as in it is determined from charmonium. For the heavy 
clover approach at (3 = 6.0 these differences are negligible within statistical errors [[IT] . 



As discussed at the beginning of this subsection, the inclusion of the terms 
(9((Aqcd/thq) 3 ) contributing to the hyperfine splitting would decrease our result by ~ 10%. 
However the heavy clover approach requires similar correction terms in the Hamiltonian to 
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achieve this level of accuracy. 

The agreement of the NRQCD and the heavy clover result for the spin-dependent hyper- 
fine splitting is encouraging, since the systematic uncertainties are quite different. The light 
quark content plays only a minor role for the hyperfine splitting, which depends essentially 
only on the heavy quark Hamiltonian. In NRQCD the leading contribution to the hyperfine 
splitting comes from the a ■ B term in the action, whereas for heavy clover this is split be- 
tween the kinetic hopping term and the clover term <J^ P T^ P , with the latter becoming more 
important as the lattice spacing becomes coarser. Both these actions give rise to systematic 
errors in the hyperfine splitting from mass-dependent radiative corrections to coefficients 
and neglected higher order terms, each at the 10% level, so the differences could have been 
significantly larger than observed. 



VII. DISCUSSION 



We present an extensive study of the B and D-meson spectrum using NRQCD heavy 
quarks and clover light quarks in the quenched approximation. 

Our results include the splitting between the strange and the non-strange meson, hyper- 
fine splittings, radially and orbitally excited states. For the first time in a lattice calculation 
we obtained a result on radially excited P-wave states. For spin-independent splittings we ob- 
serve good agreement with experimental results. However, our result for the spin- dependent 
hyperfine splitting turns out to be too low in comparison to experiment. This is a well known 
effect in quenched hadron spectroscopy. Furthermore, in the present calculation hyperfine 
splittings are also affected by the neglect of radiative corrections in the matching of lattice 
NRQCD to continuum QCD. 

Using two different values of the lattice spacing in the B-spectrum together with the 



results of [12|] allows for a detailed investigation of the residual lattice spacing dependence 
of our final results. No scaling violations are observed within the achieved accuracy. Of 
particular interest is the scaling of the B* — B s splitting, which depends heavily on the 
properties of the heavy quark content of the theory. Here scaling violations could be ruled 
out with an accuracy of ~ 10%. The P-fine structure has not been resolved for all values of 
the lattice spacings and further work is needed for this quantity. 



Our results on the B-meson spectrum are summarised in table |XXVIII| and figure 



together with the findings of |L2| . In addition to the uncertainties considered in the analysis 
at the individual values of a, the quoted uncertainties also contain an estimate for the 
residual lattice spacing artifacts of O[ot s a, a 2 ). The table gives our final results for the 
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B-meson spectrum in the quenched approximation. 

Our final results on the D-meson spectrum are shown in table |XX1X| and figure ^7 
above. This is our final result for a lattice spacing of a -1 ~ 1.1 GeV and does not include 
an estimate of the residual lattice spacing artifacts of 0{<x s a, a 2 ). For the above value of a, 
this corresponds to 13%, which is of similar size to or smaller than the otherwise achieved 
accuracy. 

We compared our results to lattice results of other collaborations obtained with NRQCD 
or in the heavy clover framework. In general we observe good agreement. Discrepancies 
which appear at first sight could be traced to underestimated errors in these other results 
or the use of different scales when converting the lattice results into physical units. The 
excellent agreement of our results with the results obtained in the heavy clover approach is 
noteworthy because of the different systematics of these approaches. 

These results are the most complete lattice results on the B and D-meson spectrum to 
date. 
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TABLE I. Simulation parameters of the gauge field configurations. For (3 = 6.2 there have been 
three different runs H, N and P with different numbers of configurations. All configurations are 
generously provided by the UKQCD collaboration. 



(3 


volume 






box size 


# configurations 


5.7 
6.2 


12 3 x 24 
24 3 x 48 




2.1 fm 
1.8 fm 


278 

H: 68; N,P: 144 


TABLE II. The hopping p 
of k c and k s are taken from [ 


ar 
17 


ameters used in 
For k s we f 


the simulation are 
live the results as 


denoted by ki to K3. The values 
determined from K, K* and cb. 


(3 K! 


K2 


K3 


K C 


K S (K) 


<s(K*) K s ((b) 


5.7 0.1380 
6.2 0.1346 


0.1390 0.1400 
0.1351 0.1353 


0.1434(1) 
0.13587(^ 


0.1399(1) 
J) 0.13466(7) 


0.1393(2) 0.1391(2) 
0.13461(+^) 0.13455(+^) 



TABLE III. Bare heavy quark masses used in the different runs at [3 = 5.7. In the second, third 
and fourth line we give the stability parameter n used in the evolution equation (||) of the runs A, 

C and S. 

qiuq 20.0 12.5 10.0 8.0 6.0 5.0 4.0 3.5 3.15 2.75 2.45 2.2 2.0 1.7 1.5 1.3 1.125 1.0 0.8 0.6 
A 1 1 1112222 3 333335 5 567 
C 1 1 1112222 3 333445 6 68 10 
S 1 2 2 3 6 8 



TABLE IV. Bare heavy quark masses and stability parameters n used in the runs H, N and P 
at |3 = 6.2. 

am Q 6.0 4.5 4.0 2.5 2.0 1.6 1.44 1.3 1.2 1.1 
HI 1 2 3 3 4 

N 1 3 3 

P 3 
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TABLE V. Smearing radii applied at (3 = 5.7 to the heavy quark Q and the light quark q. 



4>G,o 4>GJ 4>G,2 

ar Q 1.0 2.0 3.0 

ar q local local 3.0 



TABLE VI. Smearing radii applied to the heavy quarks in the run N at (3 = 6.2. The subscript 
'g' denotes a ground state hydrogenic wave function, the 'e' an excited state. Throughout this run 
we used local light quarks. 

am Q 4> Hg) l 4>Hc,l 4>Hg,2 

4.5 ar = 5.0 

2.5 ar = 5.0 ar = 8.0 

1.44 ar = 4.0 ar = 4.0 ar = 8.0 



TABLE VII. Determination of the inverse lattice spacing a -1 from the p-meson mass 
||. 17, 18, 25]. The first parenthesis gives the uncertainties arising from statistical fluctuations, the 
second the uncertainty resulting out of the chiral extrapolation. For comparison we also give the 
scales as obtained from the string tension a [26-28] and the bottomonium Xb — splitting P,|2^]. 







P 


= 5.7 




(3 = 


6.2 


Quant. 


phys. [MeV] 


lattice 


cr 1 [GeV] 


lattice 




a' 1 [GeV] 


THp 

Xb-T 


770.0(8) 
«430 
440 


0.690(8)(1° 5 ) 

0.3879(39) 

0.311(6) 


1.116(12)(lf) 

wl.10 
1.41(4)(2)(5) 


0.297(^ 2 )C 

0.1608(10) 

0.125(5) 


HO x 

-i<v 


2-59(if )(^) 

«2.67 
3.52(14) (4) (0) 
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TABLE VIII. Fitted simulation masses for the ground state pseudo-scalar and vector mesons at 
(3 = 5.7 from run A. This table has been obtained from double exponential matrix fits to correlators 
with the smearing functions 4)0,1 and 4>G,2 at source and sink. 



am.Q 




QTrisim p S 






ITUsim ,v 








k = 0.1380 


k = 0.1390 


k = 0.1400 


k = 0.1380 


k = 0.1390 


K 


= 0.1400 


20.0 


0.765(4) 


0.745(4) 


0.724(4) 


0.771(4) 


0.750(4) 


0. 


.730(4) 


12.5 


0.769(3) 


0.748(3) 


0.727(4) 


0.778(3) 


0.757(3) 


0. 


.736(4) 


10.0 


0.7715(28) 


0.750(3) 


0.728(4) 


0.7823(28) 


0.761(3) 





.740(4) 


8.0 


0.7733(26) 


0.7518(28) 


0.730(3) 


0.7868(28) 


0.765(3) 


0. 


,744(4) 


6.0 


0.7752(23) 


0.7532(27) 


0.731(3) 


0.7929(26) 


0.7712(28) 


0. 


.749(3) 


5.0 


0.7756(23) 


0.7534(26) 


0.7307(28) 


0.7966(25) 


0.7747(27) 


0. 


.753(3) 


4.0 


0.7753(22) 


0.7526(23) 


0.7296(27) 


0.8009(23) 


0.7787(26) 





.756(3) 


3.5 


0.7743(22) 


0.7514(23) 


0.7281(27) 


0.8031(23) 


0.7808(26) 


0. 


.7582(28) 


3.15 


0.7730(21) 


0.7499(23) 


0.7264(26) 


0.8046(23) 


0.7821(26) 


0. 


.7594(28) 


2.75 


0.7702(21) 


0.7468(23) 


0.7230(26) 


0.8057(23) 


0.7830(25) 





.7602(28) 


2.45 


0.7670(20) 


0.7433(22) 


0.7194(25) 


0.8062(23) 


0.7833(25) 


0. 


.7603(28) 


2.2 


0.7631(20) 


0.7392(22) 


0.7149(25) 


0.8059(23) 


0.7829(25) 





.7598(28) 


2.0 


0.7586(20) 


0.7345(22) 


0.7101(23) 


0.8050(23) 


0.7818(25) 


0. 


.7585(28) 


1.7 


0.7487(19) 


0.7242(21) 


0.6994(23) 


0.8016(23) 


0.7781(25) 


0. 


.7546(28) 


1.5 


0.7386(19) 


0.7137(21) 


0.6886(23) 


0.7969(23) 


0.7732(26) 


0. 


.749(3) 


1.3 


0.7210(19) 


0.6957(20) 


0.6702(22) 


0.7864(23) 


0.7625(26) 


0. 


.739(3) 


1.125 


0.7000(18) 


0.6743(20) 


0.6484(22) 


0.7731(23) 


0.7490(26) 


0. 


.725(3) 


1.0 


0.6788(18) 


0.6528(19) 


0.6265(21) 


0.7587(23) 


0.7343(27) 


0. 


.710(3) 


0.8 


0.6238(18) 


0.5970(19) 


0.5706(24) 


0.7183(25) 


0.6934(28) 


0. 


.668(3) 


0.6 


0.5267(20) 


0.4989(21) 


0.4709(23) 


0.6445(28) 


0.619(3) 


0. 


.593(4) 
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TABLE IX. Fitted simulation masses for the pseudo-scalar and vector meson at (3 = 6.2. The 
top section gives the results for run H, the middle one for N and in the bottom we give the result 
for run P. The results for the runs H and N have been obtained from propagators with source and 
sink smearing. In P we used source smearing with local sinks. 



am.Q 




QTTlsim.ps 






1TTLsim,v 






A 1 LAC 

K — U.I 040 


K — U.I JJ 1 


K — U.I JJJ 


A 1 1 A C 

K — U.I D4o 


K — U.I DO 1 


K — U. 1 JDJ 


o.U 


u.44y(lu ) 


U.4oo(12 ) 


u.4ly(Zl J 


U.4ou(lo ) 


u.4ol(Zu ) 


U. 410(24 ) 


4.U 


n A A o / o \ 

(J. 443(e) 


0.425(10) 


0.417(12) 


0.447(11) 


0.426(11) 


0.417(14) 


2.0 


0.420(6) 


0.405(8) 


0.398(10) 


0.430(7) 


0.412(9) 


0.404(12) 


1.3 


0.383(5) 


0.365(6) 


0.357(7) 


0.397(6) 


0.379(8) 


0.371(9) 


1.2 


0.373(5) 


0.355(5) 


0.347(7) 


0.388(5) 


0.370(8) 


0.361(9) 


1.1 


0.358(4) 


0.341(5) 


0.333(7) 


0.374(5) 


0.353(8) 


0.342(10) 


4.5 


0.435(4) 






0.440(4) 






2.5 


0.421(4) 






0.429(4) 






1.44 


0.388(3) 






0.400(3) 






1.6 


0.406(4) 






0.419(4) 
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TABLE X. Chiral extrapolation at (3 = 5.7. The first parenthesis gives the statistical uncer- 
tainty, the second one the uncertainty in the respective hopping parameter and, in the case of k c , 
the third parenthesis gives the uncertainty arising from the chiral extrapolation. 



am-Q am sim)PS am sim)V 





K c 








K C 








20.00 


0.656(6) 


(2) 


+/) 


0.726(4) (^ 7 ) 


0.663(6) 


(2) 


'-t-H\ 
■.-&) 


0.732(4) (+l 6 ) 


12.50 


0.657(6) 


(2) 


r-\-C\\ 

tl) 


0.729(4) (^ 7 ) 


0.667(6) 


(2) 


'-hU\ 

<-6> 


0.738(4) (tl 7 ) 


10.00 


0.658(5) 


(2) 


r-\-C\\ 

tj) 


0.731(4)(^ 7 ) 


0.670(5) 


(2) 


+6) 


0.742(4) (±q 7 ) 


8.00 


0.658(5) 


(2) 


+7) 


0.732(3)(^ 8 ) 


0.673(5) 


(2) 


'+V\ 

^-6) 


0.746(3) (^q 7 ) 


6.00 


0.658(5) 


(2) 


tl) 


0.733(3)(^ 8 ) 


0.678(5) 


(2) 


+5) 


0.751(3)(^o 8 ) 


5.00 


0.657(4) 


(2) 


+7) 


0.7330(28) (^ 81 ) 


0.680(5) 


(2) 


+5) 


0.755(3) (^ 8 ) 


4.00 


0.655(4) 


(2) 


+7) 


0.7319(27) (to) 


0.683(4) 


(2) 


--5) 


0.7586(29) (^ 79 ) 


3.50 


0.652(4) 


(2) 


—6> 


0.7304(26) (^o 87 ) 


0.685(4) 


(2) 


+5) 


0.7605(28) (±l so ) 


3.15 


0.650(4) 


(2) 


'+U\ 


0.7287(25) (^ 88 ) 


0.685(4) 


(2) 




0. 7617(28) (^ 81 ) 


2.75 


0.646(4) 


(2) 


'+U\ 
--6> 


0.7254(25) (^ 90 ) 


0.686(4) 


(2) 


'+U\ 
<-4> 


0.7625(28) (^ 83 ) 


2.45 


0.641(4) 


(2) 


'-t-u\ 
—6> 


0.7218(24) (+J 92 ) 


0.685(4) 


(2) 


,-4) 


0.7626(28) (^o 84 ) 


2.20 


0.636(4) 


(2) 


--6> 


0.7174(24) (t 94 ) 


0.684(4) 


(2) 


~.-4J 


0.7621(28)(^ 85 ) 


9 nn 


n fisi CK\ 

U . \JKJ ±\0 I 


(1\ 
\ z ) 


—6> 


D 71 9fif9' : iV +195N l 


n (\9,"\(a\ 

U.UOOI'i 1 


(1\ 
\ z ) 


~.-4) 


D 7fi(W98V +186N l 


1.70 


0.619(3) 


(2) 


ZD 


0.7019(23) (+<!, 98 ) 


0.678(4) 


(2) 


,-4> 


0.7569(28) (^ 89 ) 


1.50 


0.607(3) 


(2) 


ZD 


0.6911(22)(^ 01 ) 


0.672(4) 


(2) 


—4) 


0.7519(29) (^ 90 ) 


1.30 


0.587(3) 


(2) 


:i 5 °) 


0.6728(22) (^° 4 ) 


0.660(5) 


(2) 


—4) 


0.7409(29) (tl 92 ) 


1.125 


0.564(3) 


(2) 


:i 5 °) 


0.6510(21)(^ 07 ) 


0.645(5) 


(2) 


—4) 


0.727(3) 


1.00 


0.541(3) 


(2) 


;H) 


0.6291(21)(^ 10 ) 


0.630(5) 


(2) 


3) 


0.712(3)(j*>) 


0.80 


0.483(3) 


(3) 


;H) 


0.5733(23)(+^ 15 ) 


0.590(5) 


(2) 


:±S) 


0.673(4) (j*>) 


0.60 


0.380(3) 


(3) 


3) 


0.4737(22) (tf 5 ) 


0.512(6) 


(3) 


3) 


0-597(4) (j*>) 
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TABLE XI. Chiral extrapolation at (3 = 6.2. The parenthesis gives the statistical uncertainty. 
The results have been extrapolated to k c = 0.135873 and k s = 0.13466. 



crm.Q 


1 1 ^sim ,ps 




Q-TTlsim ,v 






K,. 




K,. 


D.U 


fl A 1 7/~\ C\\ 

U.41 l[L\) J 


U.44o(lU J 


U.ooo^ob J 


U.4oo(lo J 


4.0 


0.396(16) 


0.439(8) 


0.392(17) 


0.444(9) 


2.0 


0.383(12) 


0.419(6) 


0.387(17) 


0.429(6) 


1.3 


0.341(8) 


0.381(5) 


0.352(13) 


0.395(6) 


1.2 


0.331(8) 


0.371(5) 


0.342(12) 


0.384(6) 


1.1 


0.317(8) 


0.357(5) 


0.316(15) 


0.370(5) 



TABLE XII. Mass shift at (3 = 5.7. 



aniQ 


aA rc i 


aAp Crt 


20.0 


14.(5) 


18.6(11) 


12.5 


9.3(20) 


11.7(7) 


10.0 


7.9(13) 


9.4(5) 


8.0 


6.7(9) 


7.6(4) 


6.0 


5.4(6) 


5.7(3) 


5.0 


4.6(4) 


4.83(25) 


4.0 


3.96(26) 


3.89(24) 


3.5 


3.52(22) 


3.42(22) 


3.15 


3.20(19) 


3.10(23) 


2.75 


2.85(16) 


2.72(22) 


2.45 


2.57(14) 


2.44(21) 


2.2 


2.34(12) 


2.21(19) 


2.0 


2.16(11) 


2.03(17) 


1.7 


1.89(9) 


1.75(10) 


1.5 


1.68(6) 


1.56(10) 


1.3 


1.51(6) 


1.39(9) 


1.125 


1.37(5) 


1.23(7) 


1.0 


1.27(5) 


1.14(7) 


0.8 


1.16(4) 


1.01(6) 


0.6 


1.13(4) 
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TABLE XIII. Mass shift at (3 = 6.2. 



am.Q 


aA H 


^ A 

l^pcrt 


ft n 
D.U 




K QO/1 £^ 

o.oz^lo J 


A £ 

4.5 




4.4U(lzJ 


A n 

4.1) 


4.19(_ 31 ) 


o no / 1 1 \ 

3.92(11) 


2.b 




o /in/r^ 

z.4y(yj 


o n 


2 -26 (-26) 


o no /cA 


1.6 




1.64(5) 


1.44 




1.49(4) 


1.3 


1.28(1^) 


1.36(4) 


1.2 


1.17(1]]) 


1.27(4) 


1.1 


1-07(12) 


1.18(3) 
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TABLE XIV. 1-loop coefficient and q* of the perturbative expansion of the mass shift. 



am.Q 


TL 


A(1) 


ac \*A 


on nn 
zO.UU 


1 


— U.zyb8(bo J 


1 *7£. £ /^l \ 

1. Ibb{ 1 1) 


iz.oU 


1 


— U.zbUo(4b J 


1 r 7 r 7 r 7( KO\ 
L.l ( I [Ol ) 


1 n nn 
11). UU 


1 


— U.z44b(35) 


1.775(45) 


i nn 
1 .UU 


1 


— U.lyo i(oi) 


i.oU r (4y j 


£ nn 
5.UU 


1 
1 


— U.15zz(34J 


i.yby^by ) 


/i nn 
4.UU 


-i 
1 


n 1007/ 
— U.lzz /(z / J 


1.7 rfo(5b) 


A nn 
4.UU 


Z 


n 111 c^oo^ 
— U.1115(z3 J 


l.bob^bo ) 


9 £n 
o.oU 


Z 


n noon/ <~to\ 
— U.Uooy(z3 J 


1.5/4(yb) 


q nn 
O.UU 


Z 


— U.U538(zl ) 


i 97n^nn\ 

1.3 /y(yu) 


o vn 
Z. /U 


z 


n nQHQ/ on\ 
— U.UoUo(zU J 


1 07/ 1 ^ 
l.Z ) 


o ^n 
Z.oU 


z 


— U.UlU4(zo J 


U.4o(lo J 


o nn 
Z.UU 


z 


U.UOZ 1 yl6 ) 


l.zo(zo J 


i 7n 
1. (U 


z 


nil no/ 1 k\ 
U.llUy(zoJ 


1 V/1Q/'fifi\ 

l./4o(bbj 


1 fin 


z 


U.lo4b( > zoJ 


i fiQn^Q\ 
l.boy^oo J 


1.50 


2 


0.1615(24) 


1.583(42) 


1.40 


3 


0.2256(25) 


1.858(38) 


1.20 


3 


0.3276(26) 


1.765(27) 


1.00 


4 


0.5660(31) 


1.793(20) 


0.80 


5 


1.0915(42) 


1.805(13) 
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TABLE XV. Splitting between the strange and non-strange meson at (3 = 5.7. The first paren- 
thesis gives the statistical uncertainty, the second the effect of the different chiral extrapolations 
and the third the uncertainty arising from the k s determination. 



CUTlQ 


a(m pS)S 


— m.] 


») 


a(m v>s 


- m v ) 


dlTTlsav.s 


— Trig 


av) 


20.000 


0.0697(26) 


'+73\ 
v-0 > 


v-0 > 


0.0691(27) 


Uo / 


' + 1 bi\ 

Uo ; 


0.0693(27) 


f+63\ 

V-o J 


/+163\ 
v-0 > 


12.500 


0.0717(23) 


/+76\ 
v-0 > 


f+169\ 

Uo I 


0.0709(23) 


('+61 ■> 

V-0 > 


' + 1 6/ \ 

Uo 7 


0.0711(24) 


U 62n i 
V-o J 


/+'^67\ 
v-0 ) 


10.000 


0.0727(21) 


f+76\ 
v-0 ) 


(+17-\-, 

Uo ) 


0.0717(22) 


'+61 \ 

Uo , 


Uo y 


0.0720(22) 


(+ G2 \ 
V-o J 


'+170\ 
v-0 ^ 


8.000 


0.0736(20) 


r+76\ 
v-0 ) 


f+173\ 

Uo ) 


0.0725(20) 


Uo / 


/+ I'M 

Uo ) 


0.0728(20) 


f+61\ 

V-o ) 


'+172n 
v-0 ) 


6.000 


0.0749(18) 


r+74\ 
v-0 ) 


'+176\ 

Uo ) 


0.0736(19) 


(+0/ x 

Uo , 


/+ 1 /o\ 

Uo ) 


0.0739(19) 


(+58\ 
V-0 - 1 


f+~\74\ 
v-0 ) 


5.000 


0.0758(18) 


i+72\ 
v-0 ) 


f+179\ 

Uo ) 


0.0742(18) 


'+54^ 

Uo , 


Uo y 


0.0746(18) 


'+56^ 
V-0 ) 


f+'\76\ 
v-0 ^ 


4.000 


0.0770(16) 


i+7Q\ 
v-0 ) 


Uo > 


0.0751(18) 


Uo / 


, 1 1 "7^ 

'+I/A 

Uo / 


0.0756(18) 


( +53>\ 
V-0 ^ 


'+178n 

v-0 1 


3.500 


0.0779(16) 


'+69n 
v-0 > 


Uo ) 


0.0757(18) 


Uo / 


✓ I 1 75, 

Uo ) 


0.0762(17) 


'+51 \ 

V-0 ) 


'+180\ 

v-0 J 


3.150 


0.0786(15) 


r+67\ 
v-0 > 


Uo > 


0.0761(17) 


, \ Ad-, 
( "H+D^ 

Uo / 


/ 1 1 "70\ 

Uo ; 


0.0767(17) 


(+49\ 
V-0 ) 


'+181\ 

v-0 J 


2.750 


0.0795(15) 


f+65\ 
v-0 ) 


'+187^ 
Uo I 


0.0767(17) 


( +44^ 

Uo / 


s 1 1 CI X 

Uo ) 


0.0774(17) 


(+47\ 
V-o ^ 


'+182\ 

v-0 J 


2.450 


0.0803(14) 


r+63>\ 
v-0 ) 


(+189\ 

Uo ) 


0.0773(18) 


f+42^ 
Uo , 


Uo ) 


0.0780(17) 


(+45\ 
V-0 ^ 


'+184\ 

v-0 ^ 


2.200 


0.0811(14) 


v-0 ) 


(+191\ 

Uo I 


0.0778(18) 


/+4 1 •> 

Uo / 


/ + 1 OO \ 

Uo i 


0.0786(17) 


'+44n 
V-0 I 


'+185\ 

v-0 1 


2.000 


0.0818(13) 


<+59\ 
v-0 ) 


(+193\ 

Uo ) 


0.0783(18) 


/+41 

Uo / 


/+185\ 
Uo ) 


0.0792(17) 


f+43\ 

V-o 


t+\87\ 

v-0 - 1 


1.700 


0.0831(13) 


r+57\ 
v-0 1 


r+196\ 

Uo I 


0.0792(18) 


/+40^ 

Uo / 


'+187\ 

Uo i 


0.0802(17) 


f+42\ 

V-o J 


'+189\ 

v-0 ) 


1.500 


0.0842(12) 


t+56\ 
v-0 > 


f+198N 

Uo I 


0.0800(19) 


'+40^ 

Uo / 


/+188\ 

Uo ) 


0.0810(17) 


'+42n 
V-0 ^ 


(+*\9*\\ 

v-0 > 


1.300 


0.0855(12) 


f+56\ 
v-0 > 


(+202\ 
Uo ) 


0.0808(19) 


Uo , 


(+190\ 

Uo / 


0.0819(17) 


'+42n 
V-0 ^ 


r+~\93\ 

v-0 J 


1.125 


0.0869(12) 


f+56\ 
v-0 ) 


(+205\ 
Uo J 


0.0817(19) 


Uo / 


/+192\ 
Uo ) 


0.0830(17) 


'+44^ 
V-0 ) 


'+196\ 

v-0 ) 


1.000 


0.0882(11) 


'+57^ 
v-0 ) 


r+208\ 
Uo ) 


0.0824(20) 


/+4^ 

Uo , 


(+194\ 

V-0 / 


0.0839(17) 


r+45\ 
V-o ^ 


'+198\ 

v-0 J 


0.800 


0.0903(12) 


(+55\ 
v-0 ) 


'+21 3\ 

Uo J 


0.0827(22) 


/+31^ 

Uo / 


/+195x 

V-o 7 


0.0846(19) 


'+35\ 

V-o J 


'+199n 

v-0 ) 


0.600 


0.0940(12) 


f+63\ 
v-0 ) 


(+222\ 

Uo J 


0.0851(24) 


'+38^ 

Uo / 


/+201\ 

V-o i 


0.0874(20) 


f+42N 
V-0 ^ 


'+206\ 
v-0 ) 
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TABLE XVI. Splitting between the strange and non-strange meson at (3 = 6.2. Results are 
quoted for k s determined from the K meson. Fixing k s from the <b would lead to an increase by 
9%. The error bar gives the statistical uncertainty only. 



arriQ 




a(m-ps,s - m.p S ) 




a(m V)S -m v ) 


a(m saV)S - m sav ) 


6.0 




0.031(17) 




0.070(29) 


0.049(22) 


4.0 




0.043(11) 




0.052(18) 


0.050(15) 


2.0 




0.036(8) 




0.042(12) 


0.041(10) 


1.3 




0.039(5) 




0.042(10) 


0.043(9) 


1.2 




0.040(5) 




0.042(8) 


0.043(7) 


1.1 




0.040(5) 




0.055(12) 


0.051(10) 


TABLE XVII. Radially excited S- 


•wave states at (3 =5.7. These have been calculated in run S 


for k = 


0.1400. 










arriQ 




ams im (2S) 




arrisim 


(2S) - am sim (lS) 




ps 


vector 


spin-av 


ps 


vector spin-av 


8.0 


1.148(27) 


1.158(26) 


1.155(26) 


0.417(26) 


0.412(25) 0.413(25) 


4.0 


1.21(4) 


1.22(4) 


1.22(4) 


0.48(4) 


0.46(4) 0.46(3) 


3.15 


1.23(4) 


1.24(4) 


1.24(4) 


0.50(4) 


0.48(4) 0.48(4) 


2.0 


1.26(5) 


1.27(5) 


1.26(5) 


0.54(5) 


0.51(5) 0.52(5) 


1.125 


1.23(9) 


1.27(8) 


1.26(8) 


0.58(9) 


0.55(8) 0.56(8) 


0.8 


1.17(12) 


1.25(10) 


1.23(10) 


0.60(12) 


0.58(9) 0.59(9) 



TABLE XVIII. Radially excited S-wave states at (3 = 6.2. These have been calculated in run N 
for k = 0.1346. These results are extracted from double exponential vector-fits to the smeared-local 
propagators with smearing functions as listed in the second column. 



am.Q 


smearing 




amsi m (2S) 




arrisim 


(2S) - am sim i 


[IS) 






ps 


vector 


spin-av 


ps 


vector 


spin-av 


2.50 
1.44 


4>Hg,l> 4^,2 
4>Hg,1, 4>Hc,1 


0.586(32) 
0.560(32) 


0.586(36) 
0.565(34) 


0.586(35) 
0.564(32) 


0.154(34) 
0.165(33) 


0.145(39) 
0.158(35) 


0.147(37) 
0.160(34) 
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TABLE XIX. Simulation masses of the P-states at (3 = 5.7. These have been calculated in run 
S for k = 0.1400. Values given in italics are obtained from fits with low values of Q and disregarded 
in the further analysis. Therefore the spin-average has been calculated according to eq. (|l7|). The 
3 Pl and Pi operators should both yield the lightest physical J p = 1 + state. 



am.Q 


am(l 3 P ) 


am{VP-[) 


am(1 3 Pi) 


am(1 3 P 2 E) 


am(1 3 P 2 T) 


am(1P sav ) 


am(2P sav 


8.0 


1.065(28) 


1.071(26) 


1.072(27) 


1.13(8) 


1.15(7) 


1.09(7) 


1.57(21) 


4.0 


1.088(27) 


1.094(25) 


1.097(26) 


1.09(7) 


1.13(7) 


1.09(6) 


1.64(20) 


3.15 


1.094(27) 


1.103(26) 


1.105(26) 


1.08(7) 


1.12(7) 


1.09(6) 


1.68(20) 


2.0 


1.097(27) 


1.115(26) 


1.116(26) 


1.05(7) 


1.11(7) 


1.09(6) 


1.74(22) 


1.125 


1.062(28) 


1.100(26) 


1.094(26) 


1.01(8) 


1.09(7) 


1.07(6) 


1.73(24) 


0.8 


0.995(29) 


1.045(26) 


1.034(25) 


0.97(9) 


1.05(8) 


1.02(6) 


1.64(24) 



TABLE XX. Simulation masses of the P-states at (3 = 6.2. These have been calculated in run P 



for k = 0.1346 and amQ = 1.6. We report on the P-wave ground state and radially excited state. 
The spin-average has been calculated according to eq. (|l8|). 



am(1 3 Po) 




am(1 3 Pi) am(l 3 P 2 E) 


am(l 3 P 2 T) am(lP sav 


) am(2P sav ) 


0.521(17) 


0.560(19) 


0.553(20) 0.593(22) 


0.588(23) 0.568(17) 


0.90(8) 


TABLE XXI. Splittings 


of the P-states at (3 = 5.7. 


These have been calculated in run S for 


k = 0.1400. 










amQ 


am(1P-1S) sav 


am(2P-lS) sav 


am(2P-1P) sav 


am( 3 P 2 E- 3 P ) 


8.0 


0.35(7) 


0.83(22) 


0.48(19) 


0.12(10) 


4.0 


0.34(6) 


0.90(20) 


0.55(18) 


0.03(8) 


3.15 


0.34(6) 


0.93(20) 


0.59(18) 


0.00(8) 


2.0 


0.35(6) 


0.99(22) 


0.65(20) 


-0.04(8) 


1.125 


0.36(5) 


1.02(24) 


0.66(21) 


-0.07(8) 


0.8 


0.37(6) 


1.00(23) 


0.63(21) 


-0.05(9) 



TABLE XXII. Splittings of the P-states at (3 = 6.2. These have been calculated in run P for 
k = 0.1346 and am.Q = 1.6. 



am(1P-1S) sav 


am(2P-1S) sav 


am(2P-1P) sav 




am( 3 P 1 - 3 P ) am( 3 P 2 - 3 P ) 


0.152(17) 


0.49(8) 


0.33(8) 


0.039(18) 


0.032(12) 0.069(25) 
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TABLE XXIII. The hyperfine splitting rrihpf at (3 = 5.7. The directly measured results are 



obtained from the difference of the results in table VTTJ . The results for k c and k s are extracted 
from linear fits to all three simulation results. The first parenthesis gives the statistical error. At 
k c the second parenthesis gives the uncertainty of the chiral extrapolation. We quote the strange 
results for k s from the K and the second parenthesis gives the deviation of the result for k s from 
the 4>. 



CLTUq 


QiTthpf directly measured 


QiTthpf extrapolated/interpolated 




k = 0.1380 


k = 0.1390 


k = 0.1400 


K c 








20.0 


0.0054(3) 


0.0056(4) 


0.0058(4) 


0.0060(6)( 


J_1 A\ 

-o I 


0.0057(4)( 


i°) 


12.5 


0.0087(4) 


0.0089(5) 


0.0092(5) 


0.0093(8)( 


-4-91 \ 

-0 ) 


0.0090(5) ( 




10.0 


0.0108(5) 


0.0110(5) 


0.0114(6) 


0.0115(9)( 


-0 ) 


0.0111(6)( 


0) 


8.0 


0.0135(5) 


0.0137(6) 


0.0141(7) 


0.0142(10) 


Uo ) 


0.0138(7)( 


-°) 


6.0 


0.0177(6) 


0.0180(6) 


0.0185(8) 


0.0187(11) 


Uo ) 


0.0181(7)( 




5.0 


0.0210(6) 


0.0213(7) 


0.0219(8) 


0.0224(12) 


Uo ) 


0.0216(8)( 


-2 ) 


4.0 


0.0256(7) 


0.0261(8) 


0.0268(9) 


0.0276(13) 


Uo ) 


0.0265(9)( 


-\-C\ \ 

-I) 


3.5 


0.0288(8) 


0.0294(8) 


0.0301(10) 


0.0312(14) 


\-o ) 


0.0298(10) 


(-3) 


3.15 


0.0316(8) 


0.0322(9) 


0.0330(10) 


0.0343(15) 


Uo ) 


0.0327(10) 


(3) 


2.75 


0.0355(8) 


0.0363(10) 


0.0372(11) 


0.0387(16) 


/+31 \ 

Uo ) 


0.0368(11) 


(3) 


2.45 


0.0392(9) 


0.0400(10) 


0.0410(12) 


0.0427(17) 


/+31 \ 

Uo ) 


0.0406(11) 


(H) 


2.2 


0.0429(9) 


0.0437(11) 


0.0448(13) 


0.0468(18) 


U 3C h 
Uo ) 


0.0444(12) 


\-6J 


2.0 


0.0464(10) 


0.0473(11) 


0.0484(13) 


0.0506(19) 


U 3C h 
Uo ) 


0.0480(12) 


(3) 


1.7 


0.0529(11) 


0.0539(12) 


0.0552(15) 


0.0578(20) 


Uo ) 


0.0547(13) 


(3) 


1.5 


0.0584(11) 


0.0595(13) 


0.0609(16) 


0.0639(22) 


Uo ) 


0.0604(14) 


( +0 ) 
\-s) 


1.3 


0.0655(12) 


0.0668(14) 


0.0683(17) 


0.0717(23) 


Uo ) 


0.0678(16) 




1.125 


0.0731(14) 


0.0746(16) 


0.0763(19) 


0.0803(25) 


Uo ) 


0.0758(17) 


U11J 


1.0 


0.0799(15) 


0.0816(17) 


0.0833(20) 


0.0879(27) 


( +22 ) 

Uo ) 


0.0828(19) 


(+0 \ 


0.8 


0.0945(17) 


0.0965(20) 


0.0998(25) 


0.106(3)(1 


0> 


0.0990(23) 


(+0 \ 


0.6 


0.1185(24) 


0.1210(26) 


0.124(3) 


0.131(4)(± 


4 ) 

0> 


0.1230(29) 


(+0 \ 
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TABLE XXIV. The hyperfine splitting rrihpf at (3 = 6.2. The results are obtained from ratio-fits 
to the propagators with source and sink smearing. The parenthesis gives the statistical uncertainty. 


am.Q 


k = 0.1346 


k = 0.1351 


k = 0.1353 


a n 

D.U 

4.0 
2.0 
l.o 
1.2 
1.1 


n nm a( r > r >\ 

— U.UUlD^ZZ J 

0.0026(23) 
0.0077(25) 

U.UloU^Z4 J 

0.0141(24) 
0.0154(24) 


— U.UUzo^oU J 
0.0024(30) 
0.0072(32) 

0.0132(36) 
0.0146(36) 


n r\r\A (a \ 
— U.UU4(4 J 

0.0024(33) 

0.007(4) 

n m ol a\ 

U.U1Z^4 J 

0.013(4) 
0.014(4) 


4.5 
2.5 
1.44 


0.0040(5) 
0.0070(7) 
0.0120(8) 






1.6 


0.0123(18) 






TABLE XXV. Dependence of the splittings Am on the spin-averaged meson mass for (3 =5.7. 
We report the constant and linear coefficient of the dependence on 1/1X1^. For the strange split- 
tings, the number in parenthesis give the statistical error, the uncertainty from the chiral extrap- 
olation, the value of k s and the systematic uncertainty of the a-value. The last column reports 
the experimental slope from the difference of the splitting in the B and D system. Here we used 
the spin-average of the ]\ = \ states for the P-state. Uncertainties which do not apply or have not 
been evaluated for reasons detailed in the text, are marked with (-). 


Splitting 


lattice 


experiment 




Amfm^i = 0) in GeV 


S^W'-0 in GeV 2 


|^ in GeV 2 


m V)S — m v 

TTlsav,s TTLsav 


0.079(3)(i 7 )(^ 9 )(^) 
0.078(3)(^)(^ 9 )(^) 
0.079(3)(^)eM) 


0.027(8)(-)(H)(H) 
0.029(8) (-)(l 7 )(i) 


0.028(7) 

0.035(12) 

0.033(9) 


ttl P s,s(2S) - m pS)S (1S) 
m V)S (2S) -m V)S (1S) 
m saV(S (2S) -m saV)S (lS) 


0.40(3)(-)(-)(3) 
0.40(3)(-)(-)(3) 
0.40(3)(-)(-)(H) 


0.57(16)(-)(-)(H) 
0.42(16)(-)(-)(3) 
0.45(15)(-)(-)(^) 




Tn S av,s(1P) ~ m. saV)S (lS) 


0.36(8)(-)(-)(3) 


0.07(12)(-)(-)(U) 


0.11(5) 


m v — m ps 

TTlv,s TTLps,s 


-0.0001(7)(^)(-)(0) 
-0.0002(4)(-)(0)(0) 


0.151(10)(lJ 6 )(-)(iJ 6 ) 
0.144(6)(-)(i°)(lJ 5 ) 


0.297(1) 
0.304(17) 
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TABLE XXVI. Meson masses and splittings in the B system for (3 = 5.7. Overlines denote 
spin-averaged states. 

Splitting Value Uncertainties Experiment 

stat chiral strange shift a-stat a-chiral 

B s - B 85.6 MeV (20) (+£ 8 ) (+ 202 ) (3) (9) 90.2(22) MeV 

B* - B* 83.6 MeV (20) (^ 6 ) (Ij 98 ) (2) (9) (^) 91.4(38) MeV 

B s -B 84.1 MeV (20) (^ 9 ) (I 200 ) (2) (9) (± 45 ) 91.(3) MeV 

BW(2S) - BW(1S) 580.(10) MeV |,| 

B S (2S) -B S (1S) 526. MeV (38) (7) (7) (+^ 4 ) 

B*(2S) -B*(1S) 509. MeV (38) (6) (7) (+^ 2 ) 

B S (2S)-B S (1S) 513. MeV (37) - - (6) (7) (jf) - 

B*j(5850) - B S (1S) 448.(15) MeV 

B S (1P)-B S (1S) 385. MeV (70) - - (0) (4) (ff) - 

B S (2P) -B S (1P) 610. MeV (200) — — (9) (6) (^°) — 

B*-B 29.5 MeV (15) (i^ 1 ) (16) (6) (±f) 45.78(35) MeV 

B* - B s 28.3 MeV (10) — (15) (6) (jj ) 47.0(26) MeV 

B* 2 -B* 41. MeV (94) (11) (3) (ij 4 ) 
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TABLE XXVII. Meson masses and splittings in the B system for (3 = 6.2. The radially excited 



S-wave states are extracted from fits with low Q-values compared to the other results. We quote 
them in italics. 



Splitting 




Value 






Uncertainties 




Experiment 








stat 


stranj 




U oLaL 


a-chiral 




B s -B 




96. MeV 


(17) 


( + o) 




(-4) 


(12) 
v— 0/ 


90.2(22) MeV 


B* — B* 




109. MeV 


(26) 


\ — u 


) 


(-4) 


&) 

v — SJJ 


91.4(38) MeV 


B s -B 




109. MeV 


(23) 


V — u 


) 


/+3\ 
(-5) 


d) 


91.(3) MeV 


BW(2S)-B 


(*)(1S) 














580.(10) MeV 0.1 

\ / II I I II 


B S (2S) -B s ( 


IS) 


420. MeV 


(85) 




(3) 


£%) 


( + l 7 ) 




B*(2S) -B*i 


:is) 


400. MeV 


(90) 


— 


(3) 


(S) 




— 


B S (2S)-B S ( 


IS) 


405. MeV 


(90) 




(3) 


CS) 


CI 7 ) 




B*j(5850) - 


B S (1S) 














448.(15) MeV 


B s (1P)-B s l 


;is) 


395. MeV 


(45) 






{%) 






B S (2P) -B s l 


;ip) 


855. MeV 


(210) 






(il) 






B* s -B s 




27.3 MeV 


(20) 




(8) 


(3D 


Cf) 


47.0(26) MeV 


B s2~ B s0 




179. MeV 


(65) 






«) 


at) 
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TABLE XXVIII. Summary of the results on the B-meson spectrum. The table gives the average 


of our results and the result of [12 


. Errors exclude quenching 


effects but include residual lattice 


spacing artifacts of O{ot s a, or). Ag 


am overlines denote spin-averaged states. 


Splitting 


Value 


Experiment 


R R 


yu. (iu j ivie v 


on o(oo\ i\/ToV 


R* R* 


yu.(iuj iviev 


yi.4(ooj iviev 


R R 
D s — D 


yu.^iuj ivie v 


yi.^oj ivie v 


r> (lO/'TCI r(*)C1 CI 
D v ' (AO J — D J ( 1 O J 




DoU.(lUJ iviev |4Lpl] 


Rf?C) RMC1 

D [AO j — D { 1 O J 


UUU.^aUJ ivie V 




R R f 1 CI 


0411.(01)] iviev 




R*f?Cl _ R*M C) 


dzo.^ouj iviev 




R f?C1 R MCI 


530 fsn'i MpV 




R*f c ;7^?1 RMC1 




ooo.^izj iviev 


RMP) _ Rfl SI 


QUO. \ 0\J J IVlti V 




R* f^S^fil R f 1 CI 
D s jpSOUJ — D s ( \ i>) 




448.(15) MeV 


R fip) _R 1 1 CI 


411.(45) MeV 




B S (2P) -B S (1P) 


730.(200) MeV 




B*-B 


29.(5) MeV 


45.78(35) MeV 


B*-B s 


28.5(31) MeV 


47.0(26) MeV 


b 2 - B 


- 250 MeV 




B s2~ B s0 


- 250 MeV 
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TABLE XXIX. Meson masses and splittings in the D system for (3 =5.7. Overlines denote 
spin-averaged states. 



Splitting 






Value 






Uncertainties 






Experiment 










stat 


chiral 


strange 


shirt 


a-stat 


a-chiral 




D c — D 






99 9 MeV 


(13) 


(+. 2 ) 
V-0 l 


Co 34 ) 


(8) 


(14) 


V-0 / 


99 2(5) MeV 


D* — D* 






92 2 MeV 


(23") 


V-0 J 


/_|_71 Sx 

(^o ) 


(1) 


(10) 


( +50\ 
V-0 ) 


1 02 4(9) MeV 


Dc — D 






94.1 MeV 


(20) 


(+^) 


tf 1 ) 


(2) 


(11) 


V-0 1 


101 6(8) MeV 


D*(2S1 - 


n* 


1 O J 
















629 (2V6) MeV iHl 


DJ2S) - 




IS) 


665 MeV 


(130) 






(5) 


(9) 


V-0 ) 




D?(2S) - 


s 


IS) 


640. MeV 


(80) 






(9) 


(11) 


(+54\ 
V— / 




D S (2S)- 


D s 


IS) 


645. MeV 


(100) 







(8) 


(11) 


/+51\ 

V-o ) 





D**(h = 


3) - 


D s 
















483.(1) MeV 


D S (1P)- 


D s 


IS) 


411. MeV 


(61) 






(3) 


(4) 


(+27) 
V-0 J 




D S (2P) - 


D s 


IP) 


710. MeV 


(230) 






(7) 


(7) 


/+36\ 
V-0 / 




D* -D 






110. MeV 


(3) 


(3) 




(6) 


(5) 


V-0 J 


140.64(10) MeV 


D*-D s 






103. MeV 


(2) 




t 2 ) 


(6) 


(4) 


/+20\ 
V-0 J 


143.8(4) MeV 
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0.44 



0.42 



0.40 



0.38 



t — i — i — I — i — i — r 



l S ani Q = 2.5 /e = 0.1346 



¥ — i 














| e; 



w Hg | smeared— local 

□ Hg j smeared— smeared 

o Hg2 smeared— smeared 



4 



10 



FIG. 1. Dependence of the fitted pseudo-scalar simulation mass on the starting point t m i n of 
the fit range at (3 = 6.2. The results are shown for three different propagators with smearing at 
source and sink or at the source only. The octagons give those values of t m i n which, we determined 
in the Q- value analysis, to give the final result. The connecting lines are for guidance only. 
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0.00 0.05 0.10 

am q 



FIG. 2. Chiral extrapolation of the pseudo-scalar simulation mass for (3 = 5.7, am.Q = 4.0 in 
am.q = j[l — -^). Small symbols give the simulation result. The curves give a linear and quadratic 
extrapolation as described in the text. The diamond to the left gives the outcome from the linear 
extrapolation. 
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< 

cd 




6.2 



< 

td 



0.5 1 2 5 10 20 2 4 6 

am Q am Q 

FIG. 3. Comparison of the mass shift aA from the simulation and lattice perturbation theory 
as a function of the bare heavy quark mass arriQ. The octagons give the outcome of the dispersion 
relation for heavy light mesons, diamonds from heavyonium. The lines represents the perturbative 
outcome. The square gives the heavyonium shift from reference [29| for (3 = 6.2 for comparison. 
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> 

0) 



100 




in GeV~ 



FIG. 4. Mass splitting between the spin-averaged ground state S-wave heavy-light strange and 
non-strange meson. The octagons give the result for k s determined from the K-meson, squares 
from the cb-meson. For both data sets the upper errors give the uncertainty from the different 
chiral extrapolations, the lower ones the statistical uncertainty. For simplicity this figures does not 
consider uncertainties arising from the value of a, see eq. (paj) . Experimental results are displayed 
by the fancy squares. 
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B 



■b- 



- ) 



-a 



o: a = 0.076 fm 
□ : a = 0.102 fm 
O: a = 0.177 fm 
X: experiment 



0.1 



1/m. 



0.2 
in GeV 



0.3 



-l 



FIG. 5. Comparison of the strange non-strange splitting for the pseudo-scalar state at different 
values of the lattice spacing a. Octagons give our result for (3 = 5.7 and diamonds for (3 = 6.2. 



The squares give the result from [12| at (3 = 6.0 and the fancy square the experimental outcome 
for the B-meson Q. For a = 0.102 fm and 0.076 fm we give the statistical errors only, the errors 
for a = 0.177 fm have been described in figure 13. 
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1.4- am = 4.0 k = 0.1400 _I 



1.2 



l.o 



0.8 



0.6 



o 



o s 
□ m sim (2'S ) 
° m.iJl'So) 



FIG. 6. Fitting the radially excited S-state at (3 =5.7. The left hand side gives the Q-values 
and right hand side the fitted simulation masses of the pseudo-scalar ground state and the first 
radially excited state. 
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FIG. 7. Splitting between radially excited and ground state S-wave splitting at (3 
results are for k = 0.1400 w K S (K) and the spin-average. 
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FIG. 8. Fitting the 2S radial excitation at (3 = 6.2. Again we give the Q-values and on the 
right the simulation masses of the pseudo-scalar ground state and the first radially excited state. 
This example uses <bHg,i and c|)Hc,i at the sources. 
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FIG. 9. Splitting of spin-averaged P to S-wave. The results are for strange light quarks and 
the error bars give the statistical uncertainties only. The experimental result gives the B*j(5850) 
resonance and the spin-average of the D s i and D* 2 . 
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arn n = 4.0 k = 0.1400 



2.0 — 



« 1.5 



1.0 



n m S av(2P) 
O m aav (lP) 



2.0 



I 1.5 



1.0 



am„ = 1.125 k = 0.1400 



n m S av(2P) 
O m aBV (lP) 



123456 123456 

^min ^min 

FIG. 10. Fitting radially excited P-states at (3 = 5.7. The plots give the spin-averaged IP and 
IP states. The final answer is extracted from t m i n = 5 in both cases. 
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□ m sav (2P) 
° m saT (lP) 



FIG. 11. Fitting radially excited P-states at (3 = 6.2. The plots give the spin-averaged IP and 
2P state. We take our final result from t m i n = 6. 
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FIG. 12. Light quark mass dependence of the hyperfine splitting at (3 = 5.7. On the left hand 
side we show a linear fit to all 3 data points for clvtiq = 0.8. This corresponds approximately 
to the D meson. The fancy squares give the experimental result. The right hand side gives the 
fitted slope multiplied by am s , as determined from the K, for all ttlq. In order not to disguise 
the significance of our findings, the error bar gives the statistical errors of the fit parameter only. 
Experimental results are given by the fancy squares. 
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FIG. 13. The hyperfine splitting as a function of the spin-averaged heavy-light meson mass. 
Error bars give statistical uncertainties only. The fancy squares give the experimental result ||. 
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FIG. 14. The hyperfine splitting of the radially excited S-wave state is given by the octagons. 
The error bars give statistical uncertainties only. The dashed line gives the ground state hyperfine 
splitting of the strange meson for comparison. This line is not a fitted curve. 
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am Q = 2.5 /c = 0.1346 
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FIG. 15. Hyperfine splitting in the N run at (3 = 6.2. The results are extracted from the 
difference of fitted masses for the 3 Si and propagators. We display the dependence on the 
starting point t m i n of the fit range of the propagators. In all cases we used single exponential 
fits. With 'sl' we denote results obtained from smeared-local propagators, 'ss' refers to smearing 
at source and sink. The octagons give those t m j n , which we determined in the Q- value analysis to 
give the final result. The connecting lines are for guidance only. 
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0.02 



£ = 6.2 /c = 0.1346 




J I I L 



0.0 0.2 0.4 0.6 0.8 1.0 

l/am Q 

FIG. 16. Comparison of the outcome for the hyperfine splitting from the different runs at 
(3 = 6.2. The lines are for guidance only and connect the points of the different runs. 
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FIG. 17. Fit range dependence of the m( 3 P2E) — m( 3 Po) P-state fine structure from a double 
exponential matrix fit at (3 = 5.7. The upper line displays the Q-value of the fit to the 3 P2E-state. 
The fit to the 3 P -state gives Q > 0.6 even for t m | n = 1 . The bottom line gives the splitting as 
determined from the difference of the individual fit results. 
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FIG. 18. Fit range dependence of the m( 3 P2h) — m( 3 Po) P-state fine structure from the double 
exponential vector fit at (3 = 6.2. The left hand side displays the Q-value of the fit to the 3 P2E-state. 
The fit to the 3 Po-state gives even higher Q-values. The right hand side gives the splitting as 
determined from the difference of the individual fit results. 
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FIG. 19. Compilation of results for the B s — B splitting. Additional results are from p^j34 ,35|. 
Results are for k s fixed from K/p ratio, apart from the crosses. The latter use the K*/K ratio to 
fix k s , which tends to shift them downwards. Please note the bursts and crosses do not contain 
all sources of uncertainty included in the squares and fancy squares. The horizontal lines give the 
experimental result from 0. 
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FIG. 20. Scaling of the — B s splitting. The squares give our results and the fancy square the 
one from |l2] . On our finest lattice the extraction of the result turned out to be substantially more 
difficult than elsewhere, so we give this result with dashed lines, see text. The horizontal lines give 
a preliminary experimental result for an admixture of the non-strange B' — B and B*' — B* splitting 
from the DELPHI-collaboration for comparison. 
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FIG. 21. Scaling of the spin-averaged B S (1P) — B S (1S) splitting. We include the result of Q. 
The horizontal lines give the B*j(5850) resonance, which is expected to be a superposition of the 



two h 



states. 
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FIG. 22. The spin-averaged B S (2P) — B S (1P) splitting for two different values of the lattice 
spacing a. 
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FIG. 23. Results for the hyperfine splitting B* — B s for different a-values. Results from [12], 34, 35] 
are included into the plot. The bursts give results for the chirally extrapolated B* — B splitting. 
Their error bar gives only statistical errors. Crosses omit the uncertainties from the determination 
of the bare b-mass. The horizontal lines give the experimental result from |1. 
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FIG. 24. Scaling of the P-state fine structure of the B s . Results at a = 0.102 fm are from [12 
The squares and octagons are displaced for clarity. 
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FIG. 25. Results for the P-state fine structure B* 2 — B* for three different a-values. The middle 



point has been taken from [12]. Experimentally this splitting is so far unobserved. 
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FIG. 26. Spectrum of B-mesons, summarising our results and those of fl~2]| . As before, the lattice 
results, given by the octagons, give the splitting with respect to the spin-averaged IS-state B. The 
experimental results from || are included by horizontal lines. The dashed line displays a result 
from the DELPHI collaboration Q,^], interpreted to be the B^'. The P states are compared to 
the experimental result for the Bj(5732) and B*j(5850). 
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FIG. 27. Spectrum of D-mesons from our results at (3 = 5.7. Lattice results are given by 
octagons, experimental results from || by horizontal lines. The lattice results give the splitting 
with respect to the spin-average D of the IS-wave. The dashed line displays a result from the 
DELPHI collaboration Q, interpreted to be the D*'. There are two non-degenerate P-states with 
J p quantum numbers 1 + corresponding to the j\ = j an d § state. We denote these by D-] and 
respectively and similarly for the D s . The CLEO collaboration reported preliminary results for 
the Di corresponding to )\ = j @||. The shaded area gives this result. 
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FIG. 28. Comparison of recent lattice results for the strange to non-strange splitting for the 
spin-average of the D and D*-mesons. The square and both of the diamonds [38] use the K/p 
mass-ratio to fix k s , whereas the crosses |34j use the K*/K ratio. The horizontal line gives the 
experimental result from H. 
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FIG. 29. Comparison of the D s i — D s to the result of ||38[ . The experimental result for the D s i 
is again from p]. 
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FIG. 30. Comparison of our D s -meson hyperfine splitting to the findings of [34,38,3^]. Please 
note that these other results do not necessarily include all the sources of uncertainty, that we have 



included, see table XXIX . They therefore may have underestimated error bars. Again, we quote 
the experimental result from ||. See text for a discussion. 
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